Technical  Report  580^ 


M.L.  Meeks 


■) 


Radar  Propagation  at  Low  Altitudes: 

A  Review  and  Bibliography 

31  July  1981 

Prepared  for  the  Defense  Advanced  Research  Projects  Agency 
under  Electronic  Systems  Division  Contract  F19628-UO-C-OO02  by 

Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

IkXMCTON,  MASSACHVSXITS 

Approvod  for  public  relessc;  distribution  unlimited. 


81  » 


ft 


The  work  repotted  io  thii  document  wm  pttfbrniod  ar  Lincoln  Laboratory,  a 
center  for  reaeatch  operated  by  Murachuietrt  Inuitute  of  Technology.  Thii 
work  wu  sponsored  by  the  Deforce  Advanced  Retearch  Projects  Agency  under 
Ait  Force  Contract  P 19628-80-C0002  (ARPA  Order  3724). 

Thu  report  may  be  reproduced  to  satisfy  need*  of  U.a.  Government  agendet. 


The  view*  and  conclusion*  contained  in  thia  document  ate  thoee  of  the  con¬ 
tractor  and  ihouid  not  be  interpreted  at  neccMarily  tepteaenting  the  official 
poiliciea,  either  expressed  or  implied,  of  the  United  State*  Government. 


/ 

Thia  technical  report  hat  been  reviewed  and  is  approved  for  publication. 
FOR  THE  COMMANDER 

Raymond  L.  Loiselle,  Lt.Col.,  USAF 
Chief,  BSD  Lincoln  Laboratory  Project  Office 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


RADAR  PROPAGATION  AT  LOW  ALTITUDES: 
A  REVIEW  AND  BIBLIOGRAPHY 


M.L  MEEKS 
Group  48 


TECHNICAL  REPORT  580 

31  JULY  1981 

Approved  for  public  releuc;  diitribution  unlimited, 


LEXINGTON  MASSACHUSETTS 

81  9  03  U  U  3 


ABSTRACT 

'A 

This  report  reviews  the  subject  of  electromagnetic  w.ve  propagation  user 
the  eerth's  surface  ss  it  relates  to  the  detection  of  low-flying  aircraft  by 
.round-based  radars.  Separate  sections  describe  how  current  knowled,.  of  the 
three  fundamental  physical  phinomena  -  refraction,  reflection,  and  diffrac¬ 
tion  -  is  applied  to  the  problems  of  low-altitude  propagation.  Simple  models 
incorporating  these  phenomena  are  discussed:  (1)  propagation  over  a  plane 
with  arbitrary  reflection  coefficient,  (2)  propagation  over  a  knife-edge  on  a 
pirn,.,  and  (3)  propagation  over  a  smooth,  spherical  earth.  Computer  programs 
for  models  (2)  and  (3)  are  given  in  the  Appendices.  He  Include  a  bibliography 
of  the  literature  on  which  this  report  is  based:  books.  Journal  articles,  mid 
technical  reports.  References  are  listed  alphabetically  by  first  author,  mid 
subject  indices  are  given  for  major  subject  headings. 
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1. 


INTRODUCTION 


This  report  is  intended  as  a  tutorial  review  of  the  physics  of  radio 
propagation  at  low  altitudes  over  various  kinds  of  terrain.  The  objective 
here  is  to  bring  together  an  account  of  what  is  known  about  the  propagation 
effects  that  determine  the  performance  of  ground-based  radars  against  aircraft 

l 

targets  flying  at  low  altitudes. 

The  propagation  of  electromagnetic  waves  has  been  studied  actively  for 
over  50  years  with  a  number  of  objectives.  These  include  telecommunication 
and  television  coverage  prediction,  microwave- link  design,  studies  of  mobile 
radar  communication,  and  radar  performance  and  design  studies.  We  have  reviewed 
the  literature  in  these  fields,  and  we  include  here  a  bibliography  referencing, 
in  separate  sections,  books,  journal  articles,  and  technical  reports  on  this 
subject.  Document  references  in  the  bibliography  are  listed  alphabetically  by 
the  names  of  the  first-listed  author;  journal  articles  and  reports  are  addi¬ 
tionally  indexed  by  subject.  References  cited  in  the  text  distinguish  between 
books,  journal  articles,  and  technical  reports.  For  books  the  authors'  names 
appear  in  capital  letters;  for  example  BORN  ft  WOLF  (1959).  For  journal 
articles  only  the  first  letters  of  the  authors'  names  are  capitalized;  for 
example,  Day  and  Trolese  (1950) .  Finally,  for  technical  reports  the  capital 
letter  R  appears  before  the  year  designation;  for  example,  Longley  and  Rice 
(R  1968) .  Papers  that  have  appeared  in  published  proceedings  of  conferences 
have  been  considered  as  journal  articles,  but  when  conference  proceedings  are 
unpublished  or  not  generally  available,  the  papers  have  been  classified  as 
technical  reports. 

The  physical  phenomena  that  govern  propagation  —  refraction,  diffraction, 
and  reflection  (multipath)  —  are  treated  in  separate  sections.  The  final  sec¬ 
tion  describes  some  simple  propagation  models  that  combine  the  fundamental 
physical  phenomena,  and  Appendices  A  and  B  describe  computer  programs  for 
these  models. 

Notwithstanding  the  work  that  has  been  done  on  VHF,  UHF,  and  microwave 
propagation,  there  remain  a  number  of  fundamental,  unanswered  questions  that 
are  crucial  to  the  complete  understanding  of  low-altitude  propagation.  These 
questions  concern  the  reflection  properties  of  slightly  rough  surfaces  at  low 
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grazing  angles,  including  the  reflection  coefficient  of  various  kinds  of  vege¬ 
tative  ground  cover  (see  Section  4)  and  problems  associated  with  multiple  dif¬ 
fraction  effects  along  a  terrain  profile  (see  Section  S) . 


2.  THE  PATTERN  PROPAGATION  FACTOR 

In  this  section  we  formulate  the  propagation  problem  associated  with  a 
surface  radar  searching  for  incoming  aircraft  at  very  low  altitudes.  We  con¬ 
sider  a  radio  frequency  range  from  about  100  MHz  to  10,000  MHz  (VHF  through  X- 
band)  and  assume  a  geometry  in  which  the  radar  site  has  been  selected  for 
optimum  coverage  and  the  radar  antenna  mounted  at  a  height  of  30  m  or  less 
above  local  terrain.  The  aircraft,  we  assume,  will  fly  at  some  altitude  be¬ 
tween  30  and  150  m  above  ground  level.  This  search  geometry  clearly  represents 
an  extreme  case  of  low-angle  propagation. 

Propagation  effects  are  taken  into  account  in  the  radar  equation  by  in¬ 
troducing  the  pattern  propagation  factor  F.  We  designate  the  transmitted  and 
received  powers  by  Pt  and  Pr  respectively,  the  effective  aperture  and  gain  of 
the  radar  antenna  by  A  and  G,  the  range  by  R,  and  the  backscatter  cross  section 
of  the  target  by  o,  so  the  radar  equation  is 


P 


r 


Pt  AG 
(4tt  ) 2 


(2.1) 


Here  the  pattern  propagation  factor  is  defined  as 


F  ■  Electric  Field  Propagated  to  the  Target  21 

Electric  Field  at  Same  Range  in  Free'  Space  1  ‘  ‘ 

Equation  (2.1)  is  in  the  form  given  by  KERR  (1951),  page  31.  The  principle 
of  reciprocity  makes  it  unnecessary  to  distinguish  between  the  pattern  propa¬ 
gation  factor  from  radar  to  target  and  the  factor  from  target  back  to  radar. 

These  two  factors  are  equal  and  both  are  represented  by  F.  In  the  case  of  one 

2 

way  propagation  over  a  range  R,  the  received  power  is  proportional  to  F  and 
-2  4-4 

R  rather  than  F  and  R  in  the  two-way  (radar)  case.  Reciprocity  thus 
justifies  applying  the  results  of  telecommunication  studies  to  the  problem  of 
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radar  detection,  but  we  must  keep  in  mind  that,  for  example,  a  20-dB  excess 
loss,  although  acceptable  for  some  communication  purposes,  becomes  an  excessive 
40-dB  loss  in  the  radar  case. 


Referring  to  the  radar  equation,  we  can  see  that  if  propagation  over 
terrain  is  compared  to  propagation  in  free  space,  the  difference  in  return 
power  from  the  target,  measured  in  decibels,  is  given  by  40  log  F.  Alterna¬ 
tively,  if  a  minimum  detectable  return  power  Pr(min)  is  specified,  then  the 
maximum  detection  range  in  free  space,  Rmax(free  space),  is  given  by 


rt  AU 

Rmax(free  sPaco)  "  \~^r 


a 

P-(min) 


C2.3) 


\  f 

and  for  propagation  over  terrain,  the  maximum  range  becomes 


"max  ■  F  Wfree  5pa<:“)  C2'4) 

Hence,  F  is  the  factor  by  which  the  maximum  detection  range  in  free  space  is 
altered  by  propagation  effects. 

3.  REFRACTION  BFFECTS 

The  index  of  refraction  n  of  the  atmosphere  depends  on  the  pressure  P,  the 
temperature  T,  and  the  partial  pressure  of  water  vapor  e.  Over  the  frequency 
range  with  which  we  are  concerned,  the  index  of  refraction  is  essentially 
independent  of  frequency  and  given  by 

n  =  1  +  N  10"6  (3.1a) 

where 

N  -  77.6  (P/T)  +  3.73  x  105  (e/T2)  (3.1b) 


Here  P  and  e  are  expressed  in  millibars  and  T  in  degrees  Kelvin.  The  quantity 
N,  the  radio  refractivity,  is  commonly  referred  to  as  being  expressed  in  N 
units  as  given  by  Eq.  (3.1b).  The  quantities  P  and  T  can  be  measured  directly, 
but  the  water-vapor  content  is  usually  measured  indirectly  by  hygrometers 
(relative  humidity),  psychrometers  (wet-bulb  temperature),  or  dew-point  devices 
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(daw-point  temperature).  With  any  of  these  measurements,  the  saturation 
vapor  pressure  is  required  to  convert  the  measured  quantity  to  vapor  pressure. 
Figure  3.1  shows  the  saturated  vapor  pressure  of  water  plotted  as  a  function  of 
(T-273),  the  temperature  on  the  Celsius  scale.  The  contribution  of  the  water- 
vapor  component  to  the  radio  refractivity  must  be  relatively  small  for  cold  air 
because  the  saturated  vapor  pressure  is  small  as  Fig.  3.1  shows.  But  for  warm 
air  the  relative  humidity  has  a  strong  influence  on  the  value  of  N. 

Fig.  3.2  shows  the  radio  refractivity  at  sea  level  as  a  function  of  tempera¬ 
ture  for  relative-humidity  values  0,  30,  70,  and  100  percent.  We  see  that 
the  water-vapor  component  of  the  refractivity  increases  exponentially  with 
increasing  temperature. 

There  is  no  corresponding  humidity  component  in  the  atmospheric  index  of 
refraction  at  optical  wavelengths*.  For  visible  light  the  water-vapor  com¬ 
ponent  makes  a  comparatively  negligible  contribution  to  the  refractivity,  so  our 
experience  with  the  atmospheric  bending  of  light  rays  generally  does  not  apply 
in  the  radio  domain. 

Now  the  deviations  from  rectilinear  propagation  will  be  determined  by  the 
variations  in  radio  refractivity  in  the  atmosphere  through  which  the  wave 
travels.  The  dominant  variations  will  occur  in  the  vertical  direction  with  the 
pressure,  and  usually  the  temperature  and  vapor  pressure  of  water,  decreasing 
with  increasing  height  above  the  ground.  If  the  atmosphere  is  well  mixed,  that 
is  to  say  the  air  mass  has  been  thoroughly  mixed  by  convection,  eddy  turbulence, 
and  molecular  diffusion  and  is  in  a  state  of  mechanical  equilibrium  (gravita¬ 
tional  and  buoyant  forces  balanced),  then  we  can  conveniently  specify  the 
variations  in  P,  T,  and  e  with  height.  Let  us  describe  the  water-vapor  content 
in  a  well-mixed  atmosphere  by  the  specific  humidity  a,  where 

Mw 

a"  1TF  £3.2) 

w 

•The  electric  dipole  moments  of  water  molecules  can  be  reoriented  by  a  radio 
frequency  electrical  field,  but  at  optical  frequencies  the  field  reverses 
direction  too  rapidly  for  the  water  molecules  to  follow. 
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TR-580 (3.1) 


Fig.  3.1  The  saturated  vapor  pressure  of  H-0  in  millibars  vs.  temperature 
(Celsius).  Over  the  temperature  range  from  -40°C  (-40#P)  to  +40#C  (104#F)  the 
pressure  of  saturated  water  vapor  increases  by  nearly  three  orders  of  magnitude 


TEMPERATURE  CC) 


Pig.  3.2  The  radio  refractlvity  of  air  at  a  pressure  of  one  atmosphere  for  four 
values  of  relative  humidity.  As  the  temperature  rises  the  humidity  can  play  an 
increasingly  important  role  in  determining  the  refractlvity  of  the  atmosphere. 


6 


1 


\ 

\ 


} 


and  Mw  and  M  are  respectively  the  mass  of  water  vapor  and  the  mass  of  dry  air 
in  a  given  volume.  Fig.  3.3  shows  a  family  of  curves  that  specify  the  re- 
fractivity  as  a  function  of  elevation  for  several  values  of  specific  humidity 
within  the  range  frequently  encountered  in  practice. 

Fig.  3.3  reveals  three  significant  facts  concerning  the  variation  of  N 
with  elevation  in  a  well-mixed  atmosphere: 

(1)  The  slope  dN/dh,  known  as  the  refractivity  gradient  has  nearly 
the  same  numerical  value,  -0.0082  N/ft  or  -27  N/km,  throughout  the 
range  of  values  of  H  and  a  in  the  figure.  This  property  pro¬ 
vides  the  basis  for  a  simple  and  convenient  method  for  trac¬ 
ing  the  path  of  a  wave  front  through  a  well-mixed  atmosphere. 

(2)  The  curves  of  refractivity  vs.  height  are  not  quite  straight 
lines,  but  dN/dh  become  slightly  less  negative  as  height  in¬ 
creases.  However,  this  curvature  is  not  large  enough  to  affect 
the  validity  of  the  standard  ray-tracing  method. 

. (3)  The  radio  refractivity  N  is  a  strong  function  of  a  and  increases 
rapidly  as  the  specific  humidity  increases  (see  also  Fig.  3.2). 

Hence,  we  can  expect  large  changes  in  refractive  effects  in  cases 
where  the  atmosphere  ceases  to  be  well  mixed,  and  the  temperature 
is  high  enough  for  the  air  to  contain  an  appreciable  amount 
of  water  vapor  (see  Fig.  3,1). 

In  all  cases  there  is  an  upper  limit  to  the  height  of  well-mixed  regions. 

The  upper  limit  occurs  when  the  tomperature  reaches  the  saturation  value  for 
the  particular  value  of  the  specific  humidity.  This  upper  bound  can  be  recog¬ 
nized  visually  as  the  height  at  which  clouds  appear. 

Since  the  refractivity  N  decreases  with  height,  as  Fig.  3.3  shows,  we 
expect  qualitatively  that  the  effect  of  refraction  will  be  to  bend  horizontal 
rays  downward  and  carry  the  radio  waves  to  some  extent  around  the  curved  earth. 
Rather  than  repeat  the  detailed  analysis  of  the  ray-tracing  problem,  which  must 
be  solved  in  this  case,  we  shall  state  the  results  of  this  analysis.  [The 
reader  is  referred  to  LIVINGSTON  (1970),  Chapter  4,  for  a  clear  exposition.] 


C38-105 


Flg  J.3  Variation  of  rofMctlvUy  «Wh 

for  several  values  of  sp.ci  “2^1  liroujSout  th.  rang,  of  h  and  o  r.p 
tivity  dN/dh  is  «“r'.V?r.n  Iron  LIVINGSTON  (1970),  pag.  «■ 


It  turns  out  that  the  refractivity  gradient  dN/dh  causes  the  nearly  horizontal 
rays  to  be  bent  with  a  radius  of  curvature  p  given  by 


P  * 


n _ 

•  a£  008  * 


(3.3) 


where  is  the  angle  that  the  ray  makes  with  the  horizontal.  It  was  seen  in 
Fig.  3.3  that  dN/dh  is  nearly  constant  in  a  well-mixed  atmosphere,  and  so  from 
Eq.  (3.1a)  the  quantity  dN/dh  must  therefore  be  effectively  constant.  It 
follows  then  that  the  radius  of  curvature  p  must  be  constant  for  all  rays  close 
enough  to  the  horizontal  that  cos  i|i«  1.  Hence,  the  rays  in  which  we  are  in¬ 
terested  will  be  arcs  of  a  circle  of. constant  radius  p.  This  circumstance 
allows  us  to  make  use  of  an  ingenious  geometrical  transformation:  the  actual 
radius  of  the  earth  R.  (6370  km)  is  replaced  by  an  effective  radius  KR  such 
that  we  can  represent  the  refracted  rays  as  straight  lines.  The  factor  K  which 
accomplishes  this  transformation  is  given  by 


K  -  .  (3.4) 

Now  in  terms  of  the  vertical  gradient  of  index  of  refraction  n  (or  radio 
refractivity  N),  we  can  express  K  as 

K  ■  (1  *  IS1'1  ■  »  *  “>'6  »,  a?)'1  (»•« 

Here  we  have  made  use  of  the  fact  that  n  is  very  nearly  equal  to  one  as  can  be 
seen  from  Fig.  3.2.  For  a  well-mixed  atmosphere,  the  refractivity  gradient  is 
about  -27  N/km  and  K  ■  1.2.  It  has  been  customary,  however,  to  adopt  K  ■  4/3 
as  a  standard  working  value,  giving  rise  to  the  term  4/3  -  earth  atmosphere. 
But  the  well-mixed  atmosphere  would  be  a  6/5  -  earth  atmosphere.  In  practice, 
however,  radiosonde  measurements  of  atmospheric  structure  show  that  K  can 
assume  a  variety  of  numerical  values.  The  appropriate  value  of  K  will  depend 
on  geographical  location,  and  with  changing  weather  conditions  at  a  given  lo¬ 
cation,  one  can  expect  appreciable  day-to-day  changes  in  K  as  well. 
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Of  course,  atmospheric  structure  is  by  no  means  limited  to  cases  in  which 
the  refractivity  gradient  is  constant.  Temperature  inversions  and  humidity 
lapses  may  occur  to  produce  highly  variable  refractivity  profiles.  For  ex¬ 
ample,  evaporation  from  the  surface  of  the  sea  at  low  latitudes  may  lead  to 
comparatively  large  values  of  the  specific  humidity  in  a  shallow  layer  near  the 
sea  surface.  Referring  to  Fig.  3,3  we  can  understand  how  this  mechanism  could 
produce  large  values  of  N  near  the  surface  with  a  drier  otherwise  well-mixed 
atmosphere  above.  If  in  this  example  the  magnitude  of  the  refractivity  gradi¬ 
ent  becomes  large  enough,  Eq.  (3.3)  shows  that  the  radius  of  curvature  p  may 
become  sufficiently  small  for  a  ray  to  be  bent  around  the  curve  of  tho  earth  or 
even  down  to  reflect  forward  from  the  surface  of'  the  sea.  This  process  occurs 
to  produce  a  so-called  evaporation  duct  in  which  radio  waves  can  propagate 
beyond  the  horizon,  McCue  (1978)  has  given  a  complete  account  of  this  phe¬ 
nomenon  and  its  implications  for  low-altitude  propagation  at  Kwajalein. 

Many  other  atmospheric  processes  can  lead  to  ducting  and  greatly  extended 
propagation.  Basically,  ducting  can  occur  when  there  is  a  transition  from 
comparatively  high  refractivity  near  the  surface  to  distinctly  lower  refrac¬ 
tivity  values  above.  In  low  and  middle  latitudes  the  air  temperature,  particu¬ 
larly  in  summer,  may  be  high  enough  for  the  water-vapor  component  to  play  an 
important  role  (see  Fig.  3.2).  However,  when  the  temperature  is  low,  the 
atmosphere  can  contain  only  a  small  amount  of  humidity  even  when  saturated. 
Hence  in  winter  at  higher  latitudes,  ducting  must  be  a  consequence  of  tempera¬ 
ture  inversions. 

There  are  then  at  least  two  situations  in  which  extended  detection  ranges 
over  continental  land  masses  may  be  expected  due  to  refraction  effects,  Fol¬ 
lowing  rain  the  evaporation  of  water  on  the  ground  can  produce  an  over-land 
evaporation  duct.  This  phenomenon  is,  in  fact,  observed  by  weather  radars 
(private  communication,  R.K,  Crane).  Ground  clutter  is  found  to  spread  out  to 
much  greater  ranges  when  the  ground  is  wet  and  when  there  is  little  wind.  A 
totally  different  effect  may  be  expected  at  high  latitudes  on  clear  nights  over 
snow-covered  terrain.  Radiational  cooling  at  the  surface  can  produce  a  strong 
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temperature  inversion  and,  therefore,  anomalous  propagation  (Nottarp,  1967). 
This  high- latitude  phenomena  however  will  be  restricted  to  clear  nights  when 
there  is  no  wind  to  mix  the  air  near  the  ground.  But  in  both  of  these  cases, 
the  ducting  will  be  less  effective  in  increasing  detection  ranges  over  land 
than  over  water  because  the  reflectivity  of  the  ground  or  snow  surface  will  be 
less  than  that  of  a  calm  sea  or  lake  surface. 

The  widely  accepted  procedure  for  characterizing  the  effocts  of  atmos¬ 
pheric  refraction  on  radio  propagation  in  the  design  of  ground-to-ground  com¬ 
munication  systems  is  to  use  balloon-borne  radiosonde  measurements  to  estimate 
the  change  AN  in  refractivity  between  the  ground  and  an  altitude  of  100  m  above 
the  ground.  This  value  AN  for  the  first  100  meters  of  height  is  then  used  to 
estimate  dN/dh  in  Eq.  (3.5)  and,  thereby,  to  calculate  K.  The  deficiencies  in 
this  procedure  are  mainly  the  result  of  the  fact  that  a  radiosonde  is  not 
designed  to  make  accurate  low-altitude  measurements.  However,  one  can  obtain 
some  insight  into  the  variability  of  refraction  effects  at  various  geographical 
locations  because  the  radiosonde  soundings  are  usually  made  at  12-hour  inter¬ 
vals'  for  meteorological  purposes.  The  distribution  of  the  values  of  AN/Ah 
and  K  for  sites  in  Canada  and  the  northern  States  have  been  published  by  Segal 
and  Barrington  (R  1977),  and  selected  worldwide  measurements  of  refractivity 
gradients  have  been  published  by  Samson  (R  1976)  for  use  In  designing  line-of- 
sight  communication  systems. 

Pig.  3.4  shows  two  examples  of  the  probability  distributions  of  ground- 
based  refractivity  gradients  and  K  values  based  on  radiosonde  measurements  for 
Edmonton,  Alberta,  and  Inuvik,  Northwest  Territories,  taken  from  the  report  of 
Segal  and  Barrington  (R  1977).  These  plots  show  separately  the  distributions 
for  three-month  intervals  starting  in  January.  Note  that  the  values  of  K, 
reading  downward  in  these  plots,  become  increasingly  large  and  change  over  to 
negative  numbers.  This  happens  because  as  AN/Ah  becomes  increasingly  negative 
in  Eq.  (3.5)  the  expression  in  parentheses  goes  to  zero,  making  K  ■  *  when 
AN/Ah  -  -157  N/km.  For  this  value  of  the  refractivity  gradient  the  ray  curva¬ 
ture  becomes  equal  to  that  of  the  earth.  Even  more  negative  values  of  AN/Ah 
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Fig.  3.4  Two  examples  showing  the  seasonal  frequency  of  occurrence  of  refrac 
tivity  gradients  and  K  values.  The  refractivity  gradients  for  these  Canadian 
sites  in  N-units/km  were  determined  for  the  lowest  100-m  layer  of  the  atmos¬ 
phere  from  radiosonde  measurements.  These  curves  were  taken  from  Segal  and 
Barrington  (R1977) . 
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produce  negative  values  of  K,  situations  in  which  ducting  would  be  expected  to 
occur. 

From  the  examples  in  Fig.  3.4  we  can  obtain  estimates  of  the  frequency  of 
occurrence  of  various  values  of  K.  At  both  Edmonton  and  Inuvik  the  values  of 
K  lie  between  one  and  two  during  80  percent  or  more  of  the  year.  At  Edmonton, 

K  exceeds  two  or  becomes  negative  (ducting  condition)  most  frequently  during 
July-September  when  these  conditions  occur  about  20  percent  of  the  time.  At 
Inuvik  the  same  conditions  for  extended  propagation  range  (K  >  2  or  K  <  0) 
occur  during  January-March  when  the  frequency  of  occurrence  is  about  20  per¬ 
cent.  Abnormally  shortened  ranges  (0  <  K  <  1)  are  relatively  uncommon  in 
these  examples,  being  present  5  percent  of  the  time  or  less  at  both  sites. 

'Jp  to  this  point  we  have  not  considered  the  absorption  effects  of  rain 
and  atmospheric  gases  on  the  propagation  of  radio  waves.  Measurable  rainfall 
occurs  less  than  5  percent  of  the  time  except  in  coastal  areas  (see  Samson, 

R  1976),  and  rain  attenuation  has  not  been  included  within  the  scope  of  this 
report . 

A  convenient  summary  of  the  phenomena  of  rain  attenuation  has  been  given 
by  R.K.  Crane  in  Section  2.4  of  MBBKS  (1976).  Figure  3.5  shows  plots  of  the 
attenuation  one-way  in  dB/km  as  a  function  of  frequency  for  rain  (three  dif¬ 
ferent  rain  rates),  for  typical  clouds,  and  for  fog',  this  figure  appears  as 
Figure  1  (p.  178)  in  r.he  reference  cited  to  MEEKS  (1976).  For  radio  frequencies 
VHF  through  X-band  the  droplets  are  small  compared  to  wavelength  X,  and  the 
scattered  power  is  proportional  to  X’*  (Rayleigh  scattering)  when  the  droplet 
circumference  becomes  comparable  with  or  larger  than  X  the  scattering  becomes 
proportional  to  the  geometrical  cross  section  of  the  droplet  and  is  nearly 
independent  of  X  as  we  see  in  Fig.  3.5. 

Absorption  by  atmospheric  gases  is  negligible  at  frequencies  below  about 
10  GHz  compared  to  the  uncertainties  of  low-angle  propagation  even  in  the  most 
predictable  situations.  However,  significant  amounts  of  absorption  may  occur 
at  frequencies  above  10  GHz,  Absorption  by  atmospheric  gases  in  the  frequency 
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Fig.  3.5  Attenuation  coefficient  vs.  frequency  for  rain,  clouds,  and  fog 


range  from  10  to  300  GHz  is  thoroughly  discussed  by  J.W,  Waters  in  Section  2.3 
of  MEEKS  (1976). 

4.  REFLECTION  AND  ABSORPTION  EFFECTS  OF  TERRAIN 

4.1  Introduction 

In  studying  the  propagation  of  radio  waves  at  low* elevation  angles  over 
the  surface  of  the  earth  we  must  take  into  account  forward  scattering  (or 
multipath  effects).  We  consider  first  reflections  from  idealized  dielectric 
surfaces  in  the  form  of  smooth  planes  (Section  4.2).  Such  surfaces  produce 
only  coherent  reflections.  Then  we  consider  the  more  complicated  problem  of 
reflections  from  rough  surfaces  where  both  coherent  and  diffuse  scattering  can 
occur  (Section  4.3).  Finally,  we  discuss  briefly  absorption  effects  of  trees 
and  other  vegetation  which  may  cover  terrain. 

4.2  Reflection  Coefficient  for  a  Smooth  Plane  Surface 

'  1  "  ▼  ”  '  -  1  — 

When  a  wave  encounters  a  smooth  plane  surface  which  is  very  large  in 
extent  compared  to  a  wavelength,  specular  reflection  takes  place  and  the  angles 
of  incidence  and  reflection  are  equal  (Snell's  law).  For  our  propagation 
geometry  the  waves  move  in  a  nearly  horizontal  direction,  and  we  shall  use  the 
common  terminology  as  follows*,  vertical  polarization  refers  to  linear  polari¬ 
zation  with  the  electric  field  vector  lying  in  the  vertical  plane  containing 
the  incident  and  reflected  rays,  and  horizontal  polarization  refers  to  linear 
polarization  with  the  electric  field  horizontal. 

The  reflection  coefficient  T  is  defined  as  the  ratio  of  the  amplitude  of 
the  reflected  wave  to  the  amplitude  of  the  incident  wave.  The  classical  formulas 
derived  first  by  Fresnel  in  1816  give  the  reflection  coefficients  for  horizon¬ 
tal  and  vertical  polarization  as  a  function  of  the  dielectric  constant  and 
conductivity  of  the  ground  substance,  It  is  convenient  to  define  a  relative 
complex  dielectric  constant  er  as  follows: 

Er  ■  -  60  i  \a  (4-1) 
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where  e  is  the  dielectric  constant,  a  the  conductivity  in  mhos  per  meter  of  the 
ground,  eQ  the  dielectric  constant  of  free  space,  and  \  the  wavelength  in  meters. 
The  reflection  coefficient  will  furthermore  be  a  complex  quantity 

T  -  pe"1*  (4.2) 

where  p  and  $  give  the  amplitude  and  the  phase  shift  of  the  reflected  wave  with 
respect  to  the  incident  wave. 

Rather  than  repeat  the  Fresnel  equations  here,  we  shall  present  only  the 
computed  reflection  coefficients  as  a  function  of  angle  in  graphical  form  for 
representative  types  of  ground  and  for  sea  water.  Table  4.1  gives  the  appropri¬ 
ate  constants  for  several  types  of  soil  and  water.  We  have  evaluated  the  Fresnel 
equations  for  very  moist  ground,  average  ground,  polar  ice,  and  sea  water.  The 
equations  themselves  may  be  found  in  BECKMANN  AND  SPIZZICHINO  (1963),  pages 
218-219,  or  in  KERR  (1951),  pages  396-403,  The  direction  of  the  incident  wave 
is  specified  here  by  the  grazing  angle  i|i,  the  angle  between  the  incident  ray  and 
its  projection  on  the  horizontal  reflecting  plane.  Horizontal  and  vertical 
polarization  reflect  differently,  so  we  show  here  the  results  of  computations 
for  horizontal  and  for  vertical  polarization,  respectively;  the  magnitudes  of 
the  reflection  coefficients  and  py  and  their  phase  lags  <j>h  and  $y.  Figure 
4.1  shows  these  quantities  plotted  for  average  ground,  for  very  moist  ground, 
and  for  polar  ice,  which  has  nearly  the  same  parameters  as  very  dry  groiind  (see 
Table  4.1). 


Qualitatively  the  reflection  parameters  show  similar  behavior  as  a  function 
of  grazing  angle  for  various  values  of  the  surface-material  parameters.  We  will 
summarize  the  characteristics  of  the  reflected  waves  using  Figure  4.1  as  our 
example.  (1)  The  magnitude  of  the  reflection  coefficient  for  horizontal 
polarization  p^  has  the  value  one  at  very  small  grazing  angles  and  decreases 
monotonically  as  increases.  (2)  The  coefficient  for  vertical  polarization 
py  on  the  other  hand  has  a  single,  well-defined  minimum  which  is  py  ■  0  provided 
g  ■  0.  The  corresponding  grazing  angle  for  a  perfect  dielectric  is  given  by 
sin"1  (l/i/e/e0) ,  and  the  complement  of  this  angle  is  called  Brewster’s  polarizing 
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TABLE  4.1 

APPROXIMATE  ELECTROMAGNETIC  PROPERTIES  OF  SOIL  AND  WATER* 


of  Surfac< 


Seawater  at  0°C 
at  10°C 


g(mhos/m) 

4-5 

4-5 


Fresh  water  at  10°C 

84 

1  X 

10‘3  -  1 

X  10*2 

at  20°C 

80 

1  X 

10"3  -  1 

x  10‘2 

Very  moist  ground 

30 

5  x 

10“3  -  1 

X  10“2 

Average  ground 

15 

5  x 

10‘4  -  5 

x  10"3 

Arctic  land 

15 

5  x 

10'4 

Very  dry  ground  and  large  towns 

-4 

(industrial  areas) 

3 

5  x 

10'&  -  1 

x  10 

Polar  ice 

3 

2.5 

x  10"5 

*Prora  FINK Tf975) ,  pp.  18-79.  Additional  data  on  electromagnetic  properties  are 
contained  in  referendes  cited  in  the  Subject  Index  under  8.  Reflection  effects 
electrical  properties  of  surface  materials. 
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GRAZING  ANGLE  (DEG) 


Fig.  4.1  The  reflection  coefficient  as  a  function  of  grazing  angle  for  very 
moist  ground,  average  ground,  and  polar  ice.  The  relative  dielectric  constant 
and  conductivity  for  each  surface  material  are  given  in  Table  4.1.  Values  of 
p  and  <j>  in  Eq.  (4.2)  were  computed  for  frequencies  of  100  MHz  and  10  GHz, 

These  frequencies  gave  identical  curves  for  p,  but  somewhat  different  curves  for 
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angle*.  (3)  For  horizontal  polarization  the  phase  lag  $  is  nearly  tt  for  all 
values  of  i|>,  but  for  vertical  polarization  of  the  phase  lag  $  goes  from  tt  at 
small  grazing  angles  (<J»  <_  8°)  to  zero  for  large  grazing  angles  (<|i  >  45°)  with 

the  changeover  occurring  abruptly  around  Brewster's  angle.  (4)  For  very  low  grazing 
angles,  iji  in  the  neighborhood  of  one  degree  or  less,  both  and  py  are  nearly 
one  and  <j>^  and  are  nearly  it.  As  a  result,  there  should  be  little  difference 
in  the  propagation  of  horizontally  and  vertically  polarized  waves  over  ground  at 
very  low  grazing  angles.  The  conclusion  is  borne  out  by  many  propagation  experi¬ 
ments  (for  example,  LaGrone  (1960),  McPetrie  and  Ford  (1964),  and  Oxehufwud 
(1959)),  but  two  exceptions  should  be  noted:  over  a  cylindrical  ridge  of  large 
radius  of  curvature,  vertical  polarization  diffracts  somewhat  deeper  into  the 
shadow  zone  (see  Section  5,3),  and  the  trunks  of  trees  in  a  forest  scatter 
vertical  polarization  more  than  horizontal  when  the  wavelength  is  large  compared 
to  the  diameters  of  tree  trunks.  (5)  The  electromagnetic  properties  of  the 
various  ground  types  all  lead  to  effectively  the  same  reflection  coefficient  at 
very  low  grazing  angles  as  Fig.  4.1  shows. 

On  the  other  hand,  sea  water  which  has  a  comparatively  high  conductivity 
shows  somewhat  different  reflection  characteristics.  Fig.  4.2  gives  the  re¬ 
flection  coefficient  of  sea  water  at  three  frequencies  as  a  function  of  grazing 
angle.  As  one  goes  to  lower  frequencies,  the  imaginary  part  of  the  complex  di¬ 
electric  constant  of  sea  water  increases  as  Equation  4.1  specifies,  and  for  ver¬ 
tical  polarization  the  minimum  in  py  shifts  to  lower  angles  and  becomes  less 
pronounced.  Also  the  phase  lag  <f>v  makes  a  more  gradual  change  from  zero  to  it  as 
increases. 

The  special  case  of  smooth  reflecting  surfaces  discussed  in  this  section 
can  be  used  to  represent  actual  terrain  reflectivities  in  at  least  three  cases: 

(1)  ocean  or  lake  surfaces  with  negligible  wave  disturbance,  (2)  flat  desert 
surfaces,  and  (3)  snow-covered  plains.  It  turns  out  for  vertical  polarization 
that  snow  and  sand  have  quite  similar  reflection  properties  at  grazing  angles 
below  Brewster's  angles  as  shown  in  Fig.  4.3.  Nevertheless,  in  the  great 

♦This  effect  can  be  used  to  advantage  in  the  design  of  airborne  radars  which  look 
down  for  low-flying  targets.  By  employing  vertical  polarization  such  radars  can 
eliminate  or  greatly  reduce  multipath  due  to  ground  reflections. 
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Fig.  4.2  The  reflection  coefficient  as  a  function  of  grazing  angle  for  sea  ; 
water.  The  relative  dielectric  constant  is  80  and  the  conductivity  5.0  mhos/m  1 
Separate  curves  are  shown  for  frequencies  of  100  MHz,  1  GHz,  and  10  GHz.  ‘ 


C38-102 


Fig.  4.3  The  reflection  coefficient  as  a  function  of  grazing  angle  for  dry 
snow,  wet  snow,  and  sand.  The  assumed  values  for  the  relative  dielectric  con 
stant  were  1.2,  1.8,  and  3,  respectively.  A  conductivity  of  0.1  mhos/m  was 
assumed  in  all  three  cases.  Brewsters  angle  is  designated  by  \ 
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majority  of  cas  s  actually  encountered  the  ground  is  not  smooth,  and  we  must 
find  ways  to  take  into  account  not  only  the  gross  terrain  surface  features, 
but  also  the  small-scale  local  roughness  and  vegetative  cover. 

4 . 3  Reflections  from  Rough  Surfaces 

When  an  electromagnetic  wave  propagates  at  low  angles,  the  direct  wave 
combines  with  waves  reflected  from  the  ground  so  that  the  target  is  illuminated 
by  a  complex  wave  front.  We  now  consider  the  i.ature  of  ground  reflections  that 
must  be  taken  into  account  if  we  are  to  predict  propagation  under  realistic 
conditions,  that  is  over  a  rough  terrain  surface.  Let  us  think  of  the  terrain 
as  a  large  collection  of  scatterers  distributed  so  as  to  represent  the  reflec¬ 
tion  properties  of  the  surface  features.  If  we  examine  the  phase  of  the  in¬ 
dividual  waves  arriving  at  the  target  from  these  scatterers,  we  will  find,  in 
general,  that  a  fraction  of  the  phases  will  be  found  to  be  randomly  distributed 
(diffuse  scattering)  and  the  remainder  will  be  strongly  correlated  in  phase 
(coherent  scattering) .  In  the  discussions  that  follow  we  shall  neglect  the 
diffuse  scattering.  We  do  so  because  the  signals  from  the  diffuse  scatterers 
combine  with  random  phases  and  add  up  to  produce  a  field  strength  that  is  negli¬ 
gible  in  comparison  to  the  direct  wave.  On  the  other  hand,  coherent,  scattering, 
depending  on  the  nature  of  the  terrain,  can  produce  a  reflected  wave  comparable 
in  amplitude  with  the  direct  wave.  (In  fact,  if  the  terrain  should  happen  by 
chance  to  have  a  surface  of  nearly  spheroidal  contour,  the  reflected  waves  could 
be  brought  to  a  focus  at  the  target,  producing  a  reflected  wave  much  stronger 
than  the  direct  wave.)  Coherent  reflection  may  or  may  not  be  produced  depending 
on  the  roughness  of  the  terrain  and  the  wavelength  of  the  signal,  but  if  co¬ 
herent  reflection  occurs,  the  reflected  wave  will  propagate  in  the  same  azimuthal 
direction  as  the  direct  wave.  On  the  other  hand,  diffuse  reflections  will 
scatter  energy  in  all  directions.  In  the  backward  direction  the  diffuse  re¬ 
flection  of  course  produces  the  backscattered  ground  clutter  seen  by  the  radars. 

The  characteristics  of  specular  and  diffuse  reflection  are  well  described 
in  the  following  paragraphs  from  BECKMMJN  and  SPIZZICHINO  (1963),  page  241. 
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"Specular  reflection  is  a  reflection  of  the  same 
type  as  caused  by  a  smooth  surface:  it  is  directional 
and  obeys  the  laws  of  classical  optics.  Its  phase 
is  coherent  and  it  is  the  result  of  the  radiation  of 
the  points  on  the  Fresnel  elipse,  which  transmit 
waves  of  approximately  equal  phases  toward  the 
receiver.  Its  fluctuations  have  a  relatively  small 
amplitude. 

Diffuse  scattering  is  a  phenomenon  that  has  little 
directivity  and  which  consequently  takes  place  over 
a  much  larger  area  of  the  surface  than  the  first 
Fresnel  zone.  Its  phase  is  coherent  and  its  fluctuations, 
which  have  a  large  amplitude,  are  Rayleigh-distributed." 

The  existence  of  these  two  components  of  the  reflected  field  is  well  established 
by  optical  as  well  as  radio  experiments. 

A  simple  theoretical  model  has  been  developed  to  calculate  the  specular 
reflection  coefficient  R.  of  rough  terrain.  Because  of  the  assumed  randomness 

"  1  .  "  ‘  --  -  -t  .  5 

of  the  terrain,  R  must  itself  be  random  with  an  RMS  value  given  by 
s 

(Vrms  ■  (ps5RMS  r  (4,3) 

Here  pg  is  the  scattering  coefficient  and  T  is  the  reflection  coefficient  for  a 
smooth-plane  earth,  For  a  Gaussian-model  rough  surface,  according  to  BECKMANN 
and  SPIZZICHINO  (1963),  page  246,  we  have 


<  |p  |2  >  - 

''s'  average 

(4.4) 

M  .  .1.?  A»l 

(4.5) 

The  model  surface  is  specified  by  the  single  parameter,  Ah  which  is  the  standard 
deviation  of  the  normal  distribution  of  heights  measured  from  a  plane  surface. 
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ijj  is  the  grazing  angle,  and  X  is  the  wavelength.  The  quantity  A4>  is  just 
the  phase  difference  between  two  rays  that  reflect  from  the  rough  surface:  one 
ray  reflecting  from  the  mean  surface  level  and  the  other  from  a  bump  of  height 
Ah  on  the  surface. 

This  simple  theoretical  model  is  in  reasonable  agreement  with  actual 
measurements  (see  BECKMANN  and  SPIZZICHINOi  1963,  pp  317,  303)  over  a  range  of 
grazing  angles  that  goes  below  one  degree.  However,  the  propagation  measure¬ 
ments  with  which  the  model  can  be  compared  permit  only  relatively  crude  estimates 
of  Ah,  and  it  does  not  appear  useful  to  consider  more  elaborate  models  until 
more  complete  and  accurate  propagation  data  with  quantitative  ground  truth  are 
obtained. 

Notwithstanding  the  limitations  of  the  Gaussian  terrain  model,  it  is  useful 
to  examine  those  combinations  of  radio  wavelength  and  terrain  type  for  which 
this  model  predicts  appreciable  coherent  reflection.  Let  us  consider  an  ideal¬ 
ized  case  of  propagation  over  a  spherical  earth  with  homogeneous  surface  rough¬ 
ness  specified  by  Ah.  Figure  4.4  shows  a  diagram  in  which  Ah  and  X  appear  as 
the  x-  and  y-axes.  Also  indicated  in  this  figure  are  the  terrain-type  and  sea- 
state  designations  that  roughly  correspond  to  various  ranges  of  values  of  Ah. 
These  terrain  roughness  designations  are  the  ones  used  by  Longley  and  Rice 
(R  1968).  We  have  fixed  values  of  P3  and  ^  alnd  plotted  on  the  diagram  the  locus 
of  (Ah,X)  values  that  would  give  rise  to  specific  values  of  scattering  coef¬ 
ficient  at  specific  grazing  angles,  These  loci  from  Eqs.  (4.4)  and  (4.5)  must 
appear  on  the  diagram  as  diagonal  straight  lines,  all  with  the  same  slope.  We 
represent  two  cases  in  Fig.  4.4.  For  both  cases  we  assume  nominal  values  for 
the  height  of  the  radar  antenna  (30  m)  and  for  the  altitude  of  the  target  above 
ground  (90  m) .  At  a  comparatively  short  range  of  10  km  over  a  spherical  earth, 
the  grazing  angle  of  the  reflected  ray  is  0.7°,  and  at  a  longer  range  of  30  km 
the  grazing  angle  (spherical  earth)  is  0.2°.  On  the  diagram  in  Fig.  4.4  we  have 
plotted  two  bands  with  the  left  edge  of  each  band  corresponding  to  p  -  0.9 
(largely  specular  reflection)  and  the  right  edgo  to  pg  ■  0.3  (largely  diffuse 
reflection).  The  reflection  coefficient  R  in  Eq.  (4.3)  is,  of  course,  deter- 
mined  by  ps  since  at  these  low  grazing  angles  F  ■  1  (see  Section  4.2).  These 
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Pig.  4.4  Terrain-wavelength  diagram  showing  the  combinations  of  Ah  and  radio 
frequency  for  which  appreciable  specular  reflection  is  predicted  by  the  Gaussian 
terrain  model.  Two  diagonal  bands  are  shown  representing  the  transition  regions 
between  specular  and  diffuse  terrain  reflection  for  assumed  target  ranges  of 
10  and  30  km.  In  both  cases  the  antenna  and  target  heights  are  30  and  90  m. 
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bands  represent  transition  regions  between  the  uppor-left  region  in  Fig.  4.4, 
where  the  terrain  reflection  is  specular  and  the  lower-right  region  where  the 
reflection  is  entirely  diffuse.  For  the  shorter  range  of  10  km,  this  model 
predicts  that  high  reflectivity  would  be  expected  only  for  low  VHF  frequencies 
over  smooth  plains  and  for  VHF  and  UHF  frequencies  over  water.  For  the  longer 
range  of  30  km,  the  highly  reflective  regions  are  extended  to  include  slightly 
rolling  plains  for  VHF  frequencies  and  smooth  plains  for  UHF  frequencies.  The 
increased  reflectivity  for  the  longer  range  is  a  consequence  of  the  factor 
(sin  i|0  in  Eq.  (4.5).  For  all  landforms  with  which  we  are  concerned,  excluding 
very  unusual  cases,  the  value  of  Ah  will  always  be  greater  than  3m.  So,  on  the 
basis  of  this  model,  we  can  rule  out  appreciable  specular  reflection  from  ter¬ 
rain  when  the  microwave  bands  are  used  (i.e.,  frequencies  greater  than  1  QHs). 
Table  4.2  summarizes  the  types  of  terrain  that  may  produce  appreciable  coherent 
reflection  in  the  various  wavelength  bands, 

In  summary  we  can  say  that  coherent-reflection  effects  will  be  most 
Important  at  VHF  and  UHF  frequencies,  where  the  reflection  coefficients  may  ap¬ 
proach  one.  On  the  other  hand,  at  X-band,  C-band,  and  S-band  we  can  expect, 
in  almost  all  cases,  low  values  of  reflection  coefficient,  Rs  <_  0.3. 

It  should  be  emphasized  that  the  Gaussian  terrain  model  considers  p  as 

5 

a  random  variable  and  predicts  the  RMS  value  of  p  so  we  should  expect  measured 

9 

values  of  the  reflection  coefficient  to  show  fluctuations  as  the  target  position 
changes  by  small  amounts.  Since  the  model  treats  reflection  properties  of  rough 
surfaces  as  statistical  quantities,  the  most  we  can  expect  is  that  the  model 
will  predict  statistical  measures  of  Rs.  This  model,  of  course,  has  serious 
limitations,  but  nevertheless  in  the  future  evon  better  models  must  of  necessity 
be  statistical  in  character.  In  principle  if  we  had  a  very  complete  quantitative 
description  of  the  terrain,  we  should  be  able  to  predict  Rg  accurately,  but 
in  practice  the  terrain  must  be  characterized  by  statistical  quantities  unless 
we  are  willing  to  work  with  very  small  terrain  samples. 

Up  to  this  point  we  have  not  considered  specifically  the  effects  of  vege¬ 
tative  cover  on  the  reflection  properties  of  terrain.  In  general  these  effects 
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TABLE  4.2 


PREDICTED  OCCURRENCE  OF  SPECULAR  REFLECTION  (Rg  >0.5) 


Radar  Band 


Type  of  Reflection  Predicted 


X-band 

(8  -  12.5  GHz) 


Appreciable  specular  reflection  only  for  Sea  States 
<  2  (R  ■  30  km).  Overland  scattering  is  entirely 
'diffuse, 


C-band 
(4  -  8  GHz) 


Appreciable  specular  reflection  for  Sea  States 
<2  (R  ■  10  km)  and  <Z  (R  ■  30  km) .  Overland 
Tcattering  is  entirely  diffuse. 


S-band 
(2  -  4  GHz) 


Appreciable  specular  reflection  for  Sea  States 
<4  (R  *  30  km).  Overland  scattering  is  effectively 
diffuse. 


L-band 
(1  -  2  GHz) 


Appreciable  specular  reflection  for  Sea  States 
£5  and  for  land  with  Ah  <  5  m  (R  ■  30  km) ;  other¬ 
wise  diffuse  reflection. 


UHF 

(300  -  1000  MHz) 


Appreciable  specular  reflection  from  Smooth 
Plains*  and  all  Sea  States  for  R  ■  30  km. 


VHF 

(30  -  300  MHz) 


Appreciable  specular  reflections  from  Smooth  Plains* 
and  Slightly  Rolling  Plains*  for  R  «  30  km. 


*The’  terrain  types  cited  here  are  indicated  in  terms  of  Ah  in  Fig.  4.4;  these 
are  the  terrain  categories  of  Longley  and  Rice  (1968). 
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are  not  well  understood,  but  there  are  two  groups  of  measurements  that  shed 
some  light  on  the  problem.  At  a  frequency  of  4  GHz  we  have  the  results  of  measure 
ments  made  by  the  American  Telephone  and  Telegraph  Co.  during  the  engineering 
of  the  transcontinental  microwave-relay  system  (Bullington,  1954).  These 
measurements  were  made  over  a  series  of  paths  (each  roughly  25  miles  long)  that 
form  links  in  a  communication  chain  connecting  New  York  and  Denver.  In  most 
cases,  the  reflection  coefficients  fell  in  the  range  from  0. 2-0.4.  The  para¬ 
graph  below  from  Bullington’s  paper  describes  the  comparison  of  ground-reflection 
data  for  vegetated  and  barren  paths. 

"The  data  on  the  New  York-Omaha  section  of  the  route  was 
separated  from  the  Omaha-Denver  section  to  show  whether  or 
not  the  difference  in  the  type  of  terrain  adds  significantly 
to  the  spread  in  the  overall  results.  At  first,  some  of  the 
low-reflection  coefficients  were  attributed  to  absorption 
in  the  numerous  trees  along  some  of  the  paths,  but  an  ap¬ 
proximately  equal  percentage  of  low  reflections  were  found 
in  the  relatively  flat,  treeless  areas  in  eastern  Colorado." 

At  longer  wavelengths,  a  different  situation  emerges.  Extensive  measurements 
have  been  made  of  television  service  fields  at  VHF  and  UHF  frequencies.  The 
standard  receiving  antenna  height  for  these  measurements  is  30  ft,  a  height 
small  compared  to  radar  antenna  heights  or  target  heights.  Nevertheless,  these 
experimental  results  must  be  considered  relevant.  The  following  paragraphs 
are  quoted  from  a  survey  paper,  LaGrone  (1960), 

"Vegetation  is  an  important  factor  in  field-measured  data  at 
some  television  frequencies....  At  low  VHF  (frequencies),  it  was 
found  to  be  negligible  (in  its  effect),  but  at  high  VHF  (fre¬ 
quencies),  trees  and  tall  grasses  were  found  to  absorb  a 
significant  amount  of  the  signal." 

and 
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"The  absorption  of  radio  waves  by  trees  and  grasses  noted  at 
higher  VHF  (frequencies)  was  found  to  be  considerably  more 
at  UHF.  Previously  it  was  noted  that  trees  thick  enough 
to  block  vision  were  essentially  opaque  at  1000  MHz." 

The  comparatively  small  effect  of  vegetation  on  VHF  propagation  is  qualitatively 
confirmed  by  the  fact  that  at  67  MH2  no  appreciable  change  in  television  cover¬ 
age  is  found  between  summer  and  winter  conditions  in  the  Boston  area  (private 
communication,  Joseph  Blake,  Chief  Engineer,  WBZ-TV,  Channel  4).  Clearly  more 
experimental  work  will  be  required  to  understand  the  effects  of  vegetative  cover, 
particularly  at  UHF  frequencies,  for  antenna  (or  target)  heights  from  100  ft  to 
500  ft,  rather  than  30  ft. 

5 .  DIFFRACTION  EFFECTS 

0. 1  Introduction 

It  is  the  nature  of  waves  to  bend  around  obstacles.  This  diffraction 
phenomenon  plays  a  central  role  in  the  propagation  of  radio  waves  at  low  ele¬ 
vation  angles  over  the  ground.  Hills  and  ridges  diffract  energy  into  the 
shadow  zone  and  make  possible  radio  communication  and  radar  target  detection 
within  these  zones,  Theoretical  models  that  describe  diffraction  by  a  knife- 
edge  and  by  cylinders  with  various  radii  of  curvature  may  be  used  in  predicting 
radio  propagation,  and  the  areas  of  application  and  limitations  of  these  simple 
geometrical  models  are  fairly  well  understood.  In  this  chapter  we  discuss 
these  models. 

In  the  theoretical  analysis  of  diffraction,  the  obstacles  are  frequently 
chosen  to  be  perfect  conductors  and  the  polarization  of  the  elec  /omagnetic 
waves  may  be  an  important  consideration.  However,  for  knife-edg 3  diffraction 
the  results  are  rigorously  independent  of  polarization,  and  for  cylinder  dif¬ 
fraction  the  polarization  dependence,  although  strong  for  perfect  conductors, 
is  small  for  the  values  of  dielectric  constant  and  conductivity  encountered  on 
terrain  at  our  frequencies  of  interest.  In  this  chapter  we  shall  assume  that 
the  polarization  is  horizontal  (i.e.,  the  electric  vector  is  perpendicular  to 
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the  plane  containing  the  incident  and  diffracted  rays).  For  diffraction  by 
cylinders,  Hacking  (1970)  has  confirmed  experimentally  that  vertical  polariza¬ 
tion  propagates  somewhat  more  strongly  into  the  shadow  zone  than  horizontal 
polarization,  but  the  difference  is  small  (<2  dB) . 

We  shall  consider  in  the  sections  that  follow  the  Fresnel -Kirchhoff  cal¬ 
culation  of  diffraction  by  a  knife-edge  and  the  corrections  that  must  be  added 
when  it  becomes  necessary  to  account  for  the  finite  radius  of  curvature  of  the 
diffracting  mask.  It  turns  out  that  most  hills  and  ridges  can  be  well  repre¬ 
sented  by  a  knife-edge,  In  the  final  section  we  discuss  the  problem  of  successive 
diffraction  over  several  knife-edges. 

The  diffraction  models  which  we  discuss  in  this  chapter  are  essentially  two- 
dimensional  models,  so  a  question  arises  concerning  the  widths  that  obstructions 
must  have  perpendicular  to  the  plane  of  incidence  for  these  models  to  apply. 

The  answer  to  this  question  is  that  the  obstruction  must  have  a  transverse 
width  greater  than  the  f irst-Fresnel-zone  width, 

[1  -  (^p)  -  1)2]1/2 

where  d  is  the  distance  from  the  radar  to  the  target  (or  transmitter  to  receiver), 
and  d^  is  the  shorter  of  the  two  distances:  radar  to  obstruction  and 
obstruction  to  target.  For  an  obstruction  centered  on  a  10-km  path,  the  first - 
Fresnel-zone  widths  are  87  m,  27  m,  and  9  m  for  radio  frequencies  100  MHz,  1  GHz, 
and  10  GHz,  respectively. 

5 . 2  Diffraction  by  a  Knife-Edge 

The  problem  of  diffraction  of  an  electromagnetic  wave  by  a  knife-edge  is  a 
classical  problem  in  optics  solved  by  Sommerfield  in  1896,  and  the  detailed 
solution  is  given  in  BORN  8  WOLF  (1959),  Sections  8.7  and  11.5.  The  geometry 
of  the  problem  is  illustrated  in  Fig.  5.1, 

Here  we  consider  first  one-way  propagation  from  a  transmitter  at  T  to  a  re¬ 
ceiver  at  R  over  a  knife-edge  at  0.  An  important  parameter  is  the  clearance 
of  the  line  TR  over  the  knife-edge  or,  if  the  line  TR  intersects  the  mask, 
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the  distance  from  0  to  TR.  We  call  this  clearance  distance  A,  taking  A  as 
positive  when  TR  clears  the  mask  and  as  negative  when  TR  intersects  the  mask 
(R  in  the  shadow  region).  Figure  5.2  shows  this  definition  of  A. 

The  solution  to  the  knife-edge  diffraction  problem  gives  the  field  at 
R  due  to  a  source  at  T,  and  from  the  principle  of  reciprocity  the  same  field 
would  be  produced  at  T  if  the  transmitter  were  placed  at  R.  Accordingly,  we 
consider  the  transmitter  and  knife-edge  as  fixed  and  denote  the  field  at  the 
receiver  as  E(dR,  yR).  The  explicit  solution  E(dR,  yR)  for  wavelength  A 
near  the  geometric  shadow  (for  either  polarization)  is  given  by 


E(dR,  yR) 


e-i(ir/4) 

rr 


ikr 


Uj  +  C(w)]  +  i  [£  +  S(w) ] } 


(5.1) 


where: 


k  -  2rr/A 

r  *  ^(yT  -VR)2  *  (dT  ♦  dR)“ 

w  -  vT"  ( - £ - ) 

w 

C(w)  =  /  cos  j  ‘c2  dr 
a 

W  - 

S(w)  =  /  sin  y  r  dt 


(5.2) 


o 

The  functions  C(w)  and  S(w)  are  the  Fresnel  integrals  and  their  argument 

w  has  a  simple  geometrical  interpretation.  When  the  clearance  A  has  the 

value  A^  such  that 
o 


A0  -  /A  d^,  d r7  (drp  '+'  <iR) 


(5.3) 


then  the  path  length  TOR  is  longer  than  the  direct  path  TR  by  an  amount  A/2. 
Thus,  A0  as  defined  above  is  just  the  clearance  that  puts  the  knife-edge  at  the 
boundary  of  the  first  Fresnel  zone,  and  from  Eq.  (5.2)  the  Fresnel  parameter 
w  may  be  expressed  as 
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w  b  clearance  of  ray  TR _  A  f5  41 

yry  "  clearance  of  lirst  Fresnel  zone  ’  3^ 

In  other  words,  w/vT-  is  just  the  clearance  of  ray  TR  expressed  in  units  of 
the  first-Fresnel-zone  clearance.  (A  more  detailed  discussion  of  this  problem 
is  given  in  BRODHAGE  $  HORMUTH  (1977),  Section  3.1.) 

Now  let  us  look  at  the  solution  to  the  knife-edge  diffraction  problem  given 
by  Eq.  (5.1).  Figure  5.3  shows  the  power  propagated  from  T  to  R  relative  to 
free-space  propagation.  This  quantity  expressed  in  decibels  is  equal  to  20  log  F, 
where  F  is  the  pattern  propagation  factor  defined  in  Chapter  2. 

The  independent  variable  in  this  figure  is  just  the  ray  clearance  in  units 
of  A0.  Expressed  in  this  way  the  result  is,  of  course,  independent  of  wavelength. 

In  order  to  see  how  radio  waves  of  different  frequencies  diffract  over 
a  knife-edge,  let  us  consider  an  example.  Figure  5.4  shows  the  geometrical 
arrangement  and  a  plot  of  the  pattern  propagation  factor  F  computed  as  a  function 
of  the  target  height  for  five  different  radar  frequencies  (wavelengths):  100  MHz 
(3  m),  300  MHz  (1  m) ,  1  GHz  (30  cm),  3  GHz  (10  cm),  and  10  GHz  (3  cm).  It 
is  convenient  to  interpret  the  quantity  F  as  the  factor  that  multiplies  the 
maximum  range  of  the  radar  in  free  space  to  obtain  the  maximum  range  under 
propagation  conditions  characterized  by  F  (see  Section  2,  Eqs.  (2.3)  and  (2.4)). 
Hence,  a  radar  with  a  maximum  range  in  free  space  of  150  km  would  be  able  to 
detect  the  target  in  Fig.  5,4  at  a  range  of  15  km  provided  F  ^  0.1  in  the  ab¬ 
sence  of  ground  clutter.  The  curves  plotted  in  Fig.  5.4  show  that  the  target 
would  be  detectable  at  any  height  above  the  ground  at  radar  frequencies  below 
300  MHz.  However,  for  higher  radar  frequencies  --  1,  3,  and  10  GHz  --  the  mini¬ 
mum  heights  for  detection  are  62,  91,  and  100  m  above  the  ground,  respectively. 
This  is,  of  course,  a  simplified  example  which  neglects  ground  reflection  and 
ground  clutter.  In  Section  6  we  consider  the  combined  propagation  effects 
of  diffraction  and  reflection. 

In  the  design  of  microwave  links,  careful  consideration  must  be  given  to 
the  diffraction  effects  of  terrain.  Field  experiments  (for  example,  Oxehufwud, 
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Fig.  5.3  Propagation  over  a  knife-edge  as  a  function  of  normalized  clearance. 
The  abscissa  represents  the  clearance  A  normalized  in  units  of  first-Fresnel- 
zone  clearance.  The  ordinate,  20  log  F  (F  -  pattern  propagation  factor,  see 
Chapter  2),  is  also  the  one-way  propagated  power  relative  to  free  space  expres¬ 
sed  in  dB. 
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Pig.  5.4  An  example  of  knife-edge  diffraction  for  five  radar  frequencies.  The 
absicssa  of  the  graph  represents  target  height  as  shown  in  the  diagram.  The 
ordinate  represents  both  20  log  P  and  P,  where  F  is  the  pattern  propagation 
factor.  The  dashed  horizontal  line  shows  the  detection  limit  for  targets  in 
the  shadow  region.  Above  this  line  a  target  could  be  detected  (in  the  absence 
of  clutter)  by  a  radar  with  a  maximum  free-space  range  of  150  km. 
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1959)  have  shown  that  terrain  obstruction  will  not  be  evident  if  60  percent 
or  more  of  the  first  Fresnel  zone  is  unobstructed.  Fig.  5.3  shows  that  for 
A/Aq  0.6  the  effect  of  knife-edge  diffraction  on  one-way  propagation  would 
be  no  more  than  l.S  dB.  According  to  BROOHAGB  4  HORMUTH  (1976),  Section  4, 
the  design  practice  is  to  place  the  antennas  high  enough  that  the  first  Fresnel 
zone  is  entirely  unobstructed  in  order  to  allow  for  unfavorable  atmospheric 
refraction  conditions.  Depending  on  the  refractivity  gradient  in  the 
atmosphere,  the  earth-radius  factor  K  (see  Section  3)  will  determine  the 
effective  clearance.  In  constructing  computational  models  of  low-angle  propa¬ 
gation,  one  therefore  has  the  support  of  theory  and  practice  in  neglecting 
diffraction  for  cases  in  which  A/AQ  £  0.6.  In  such  cases  we  have  what  is  re¬ 
ferred  to  as  Fresnel  clearance. 


5.3  Diffraction  by  Cylinders 


In  some  situations  a  more  appropriate  model  for  diffraction  by  a  hill  or 
ridge  is  provided  by  the  case  of  diffraction  by  a  cylinder.  Solutions  to  the 
problem  of  diffraction  of  electromagnetic  waves  around  a  cylinder  have  been 
given  by  Rice  (1954),  Neugebauer  and  Bachynski  (1958),  and  by  Wait  and  Conda 
(1959).  We  shall  make  use  of  the  results  of  Dougherty  and  Mahoney  (1964),  who 
extended  the  solution  of  Wait  and  Conda  (1959)  and  presented  the  results 
in  convenient  graphical  form.  Cylinder  diffraction  is  presented  in  the  form 
of  corrections  to  be  applied  to  the  solution  for  a  knife-edge.  Dougherty 
and  Mahoney  (1964)  use  the  dimensionless  parameter  p  to  characterize  the 
effect  of  the  finite  radius  of  curvature  of  the  cylinder.  They  define  p  as 

.VTT" 

p  -  V  —z —  V  -j— T— 


(S.5) 


where  r  is  the  radius  of  curvature  of  the  cylinder  and  A  is  the  wavelength. 

Fig.  5.5  shows  the  diffraction  loss  as  a  function  of  the  normalized  clearance 
A/Aq  for  six  different  values  of  p.  The  curve  for  p  ■  0  corresponds  to  the  pre¬ 
vious  case  of  knife-edge  diffraction.  For  a  frequency  of  300  MHz,  for  example, 
and  a  hill  with  radius  of  curvature  0,5  km  forming  an  obstruction  midway  along 
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Fig,  5.5  Propagation  over  a  cylinder  as  a  function  of  normalised  clearance. 
The  parameter  p,  defined  in  Eq.  (5.5),  increases  with  increasing  radius  of 
curvature  of  the  cylinder.  Knife-edge  diffraction  corresponds  to  p  •  0. 
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a  10-km  path,  we  calculate  from  Eq.  (5.5)  the  value  p  =  0.131.  Referring  to 
Fig.  5.4,  we  can  see  that  for  h/hQ  ■  -0.5  the  one-way  diffraction  loss  is 
12  dB  for  a  knife-edge  and  14  dB  for  the  rounded  hill.  So,  in  this  case,  the 
rounded  hill  gives  2  dB  more  loss,  a  comparatively  small  change.  For  shorter 
ranges,  Eq.  (5.5)  would  give  a  larger  value  of  p  and  greater  excess  loss  for 
an  obstruction  of  this  shape.  If  we  go  to  higher  frequencies  (shorter  wavelength) 
we  would  find  little  decrease  in  p  because  X  enters  only  to  the  power  1/6;  also 
the  dependence  on  r  goes  as  the  power  1/3,  which  means  that  p  is  also  not  a 
strong  function  of  r. 

However,  when  the  radius  of  curvature  begins  to  approach  an  earth  radius 
with  the  corresponding  p  >  1,  the  excess  loss  over  that  of  a  knife-edge  may  be¬ 
come  very  large.  This  circumstance  leads  to  the  concept  of  obstacle  gain  where¬ 
by  propagation  may  actually  be  improved  by  erecting  a  knife-edge  barrier  near 
the  center  of  a  spherical -earth  path,  For  example,  Dickson,  et  al.  (1953), 
reported  that  for  one-way  propagation  on  a  160-mile  path  over  the  8,000-ft 
level  of  Mt,  Fairweather  in  Alaska  at  a  frequency  of  38  MHz,  the  measured  loss 
was  73  dB  less  than  that  calculated  from  smooth-earth  diffraction  theory,  an 
"obstacle  gain"  of  73  dB,  At  VHF  frequencies,  according  to  Dickson,  et  al.  (1953), 
it  is  commonly  observed  that  very  long  paths  over  mountainous  terrain  provide 
good  communication  with  negligible  fading. 

5 . 4  Multiple  Diffraction 

Radio  waves  propagating  near  the  surface  of  the  earth  may  frequently 
encounter  more  than  one  diffracting  mask.  The  problem  of  multiple  diffraction 
has  been  investigated  by  Millington,  Hewitt,  and  Immirzi  (1962)  in  the  case 
of  double  knife-edges;  their  analysis,  based  on  Fresnel  diffraction  theory, 
cannot  easily  be  generalized  to  more  than  two  successive  knife-edges.  However, 
Deygout  (1966)  has  developed  a  semi-empirical  procedure  for  predicting  multiple 
diffraction  losses  which  agrees  with  the  results  of  Millington  et  al.  for  the 
double  knife-edge  and  gives  good  agreement  with  propagation  measurements  over 
paths  with  up  to  five  diffracting  hills. 


The  Deygout  method  of  calculation  as  reported  in  Deygout  (1966)  gives 
significantly  better  agreement  with  experiment  than  the  methods  of  Bullington 
(1947)  and  Epstein  and  Peterson  (1953). 

We  shall  describe  the  Deygout  method  by  showing  an  example  in  which  three 
successive  diffractions  are  handled.  The  method  of  computation  can  be  under¬ 
stood  as  a  straightforward  general!  r;«Hcn  of  this  example.  Pig.  5.6(a)  shows 
the  clearances  of  the  three  masks  f.win  the  line  TR.  These  clearances  are 
by  definition  negative  numbers,  The  principal  mask  must  be  determined  by 
dividing  the  clearance  h  for  each  mask  by  its  first-Fresnel-zone  clearance 
for  path  TR  and  selecting  the  most  negative  Fresnel  clearance  as  the  principal 
mask.  In  this  example  M0  is  the  principal  mask.  Next,  we  draw  the  paths  to  the 
the  top  of  the  principal  mask,  TMj  and  M^R,  and  record  A2  ■  h2  as  the  clearance 
for  the  principal  mask.  Fig.  5.6(b)  shows  this  step  in  the  construction.  Finally, 
we  draw  paths  TM1  and  and  paths  MjMs  and  MjR  as  shown  in  Fig.  5.6(c),  and 
the  appropriate  clearances  A^  Aj,  and  Aj  are  labeled.  Note  that  Aj  is  now  a 
positive  clearance.  Now  the  excess  path  attenuation,  one-way,  can  be  computed 
by  taking  into  account  the  clearances;  A1  for  path  TM2,  A2  for  path  TR,  and 
Ag  for  path  M2R.  The  corresponding  losses  in  decibels  can  be  determined  from 
Fig,  5.3  and  added  together  to  obtain  the  total  excess  loss, 

Errors  in  the  Deygout  method  of  calculation  may  be  estimated  by  comparing 
the  attenuations  calculated  for  two  successive  knife-edges  by  the  Deygout 
method  and  the  more  rigorous  method  of  Millington,  ot  al.  (1962).  This  compari¬ 
son  shows  that  the  Deygout  method  overestimates  the  attenuation  when  the  two 
knife-edges  have  nearly  the  same  Fresnel- zone  clearance  considered  separately. 

More  specifically,  errors  will  occur  when  the  secondary  mask  has  a  Fresnel- zone 
clearance  that  is  70  percent  or  more  of  the  Fresnel-zone  clearance  of  the 
principal  mask.  If  the  two  nearly  equal  masks  are  close  together,  the  over¬ 
estimate  produced  by  the  Deygout  method  may  approach  6  dB.  Therefore,  one  must 
be  careful  when  employing  the  Deygout  method  in  cases  where  two  or  more  promi¬ 
nent  features  in  the  terrain  are  close  together  and  have  nearly  the  same  heights. 


from  which  the  principal  mask  if  determined h  #  *how*  the  clear*nce 
to  the  smaller  masks,  and  (J)  ahowthl^iiviJ3  ahow*  the  P«ths  with  reap 
Attenuations  of  the  three  masks  ™  clwr*nce,  u*®<*  to  calculate  the  P 
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PROPAGATION  MODELS 


6.1  Introduction 

Hero  we  consider  the  problem  of  computing  the  signal  strength  propagated 
over  a  specific  terrain  profile  taking  into  account  the  effects  of  refraction, 
reflection,  and  diffraction.  Following  the  discussion  of  refraction  in 
Section  3,  we  see  that  in  most  cases  the  downward  bending  of  ray  paths  by  the 
atmosphere  can  be  taken  into  account  by  substituting  for  the  actual  curvature 
of  the  earth  a  modified  radius  of  curvature  that  is  somewhat  larger.  The 
factor  4/3  is  commonly  considered  as  the  standard  correction  factor,  but  over 
land  this  factor  varies  with  weather  conditions  lying  between  one  and  two 
roughly  90  percent  of  the  time  (see  Section  3).  Since  maps  describe  the 
terrain  relief  with  reference  to  sea  level,  one  can  correct  the  sea  level 
elevations  for  earth  curvature  and  refraction  by  adding  a  term  as  follows: 


i 


hi  "  hi(map)  ”  TTT 


where:  h 


h,  (map) 


.1 

K  R 


the  height  of  an  element  of  the  terrain  profile, 
the  map  height  above  sea  level, 
the  distance  from  the  radar  to  the  element, 
the  modified  earth  radius. 


In  this  way  the  terrain  profile  may  be  obtained  with  the  refraction  correction 
built  in. 


The  influence  of  terrain  reflection  is  not  so  easily  taken  into  account. 
The  discussion  in  Section  4  points  out  that  appreciable  coherent  reflection 
from  the  ground  can  be  expected  only  for  terrain  classified  as  plains,  and  pri¬ 
marily  for  radars  operating  in  VHP  or  UHF  frequency  bands.  However,  the 
magnitude  of  the  reflection  coefficient  will  depend  on  the  vegetative  ground 
cover,  and  we  lack  an  experimental  basis  for  assigning  reflection  coefficients 
to  common  types  of  ground  cover,  Nevertheless,  we  can  examine  the  effects  of 
ground  reflection  in  some  particular  cases  by  trying  several  values  for  the 


reflection  coefficient.  In  all  terrain  types  except  plains,  diffraction  by 
hills  and  ridges  will  be  the  most  important  phenomenon  governing  radio  propa¬ 
gation  at  low  altitudes.  When  there  is  a  single  diffracting  ridge  or  hill 
along  the  path,  the  problem  is  easily  handled  with  the  methods  described  in 
Section  5,  If  there  are  several  diffracting  masks  in  series  along  the  profile, 
the  propagation  may  be  calculated  by  the  Deygout  method  which  gives  a  good 
approximation  in  most  cases  (see  Section  S.4).  On  the  other  hand,  over  smooth 
plains,  radio  waves  reflect  from  and  diffract  over  what  is  essentially  a  sphere 
(or  cylinder)  with  a  very  large  radius  of  curvature,  roughly  that  of  the  earth. 
This  problem  has  been  solved  in  general  only  for  a  smooth  dielectric  or  conducting 
sphere. 

In  the  sections  that  follow  we  examine  several  specialized  propagation 
models:  (1)  propagation  over  a  plane,  (2)  propagation  over  a  knife-edge  on  a 
plane,  and  (3)  propagation  over  a  smooth  spherical  earth,  The  solutions  to 
these  problems  give  some  insight  into  the  problems  of  propagation  modeling. 

6.2  Propagation  Ovqr  a  Plane 

The  standard  method  of  investigating  propagation  over  a  flat  plane  makes 
use  of  the  fact  that  the  plane  reflects  like  a  mirror,  and  the  reflected  wave 
can  be  assumed  to  originate  at  the  mirror  image  of  the  source.  In  effect, 
one  replaces  the  reflecting  plane  by  the  image  of  the  source  and  assigns  a 
phase  shift  to  the  image  wave  to  take  into  account  the  phase  change  upon 
reflection.  Also,  the  amplitude  of  the  reflected  Wave  must  be  multiplied  by 
a  factor  equal  to  the  magnitude  of  the  reflection  coefficient.  The  direct 
wave  and  the  wave  from  the  image  combine  to  form  a  pattern  of  lobes  and 
nulls.  Figure  6.1  shows  the  geometry  of  the  problem,  If  the  horizontal 
distance  to  the  target  is  large  compared  with  the  heights  z ^  and  z2,  then  the 
path  difference  is: 

AR  ■  dj  -  d^  ■  2  z^  Zj/x 

Let  us  assume  that  the  source  is  a  radar  with  the  axis  of  its  antenna  beam 
pointed  horizontally  toward  the  target  and  that  the  pattern  of  the  antenna  is 
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given  by  f(0),  where  8  is  the  angle  between  the  horizontal  and  the  target 
viewed  from  the  radar  antenna  or  from  the  mirror  image  of  the  antenna.  Then 
the  pattern  propagation  factor  will  be  given  by: 

F  -  | f (GjD  +  pf(62)  e“i(2TT  *R/X  * 

where:  R  *  range  from  radar  to  target, 

-  (z2  -  z^/R  (radians), 

>  ©2  ■  (ij  +  z  2 ')  /  R  (radians), 

pe"^  ■  complex  reflection  coefficient, 

X  ■  wavelength. 

Since  the  phase  shift  will  normally  be  it  at  low  elevation  angles  (see 
Section  4.2),  nulls  will  be  produced  when  the  target  height  has  values: 

z2  “  °»  2(^'  3(^  1  ‘  * 
and  peaks  produced  for  target  heights 

z2  “  *  1  * 

The  widths  of  the  lobes  are  thus  directly  proportional  to  wavelength  and  in¬ 
versely  proportional  to  radar-antenna  height. 

If  we  neglect  the  antenna  pattern,  then  the  ratio  of  the  null  signal 
strength  to  the  peak  signal  strength,  expressed  in  decibels,  is  20  log  (l-p)/(l+p). 
For  example,  the  null-to-peak  power  ratios  are  -15  dB  for  p  *  0.7,  -9.5  dB 
for  p  ■  0.5,  and  -5.4  dB  for  p  ■  0.3.  A  standard  method  for  measuring  reflec¬ 
tion  coefficient  makes  use  of  the  relationship  between  null-to-peak  ratios 
and  reflection  coefficient. 

6.5  Propagation  Over  a  Knife-Edge  on  a  Plane 

The  problem  of  propagation  over  a  flat  plane  can  be  generalized  to  in¬ 
clude  a  knife-edge  mask  on  the  plane.  This  is  accomplished  by  using  images 
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to  represent  the  reflected  waves.  Since  reflections  may  occur  on  both  sides 
of  the  knife-edge,  images  of  both  the  radar  and  the  target  must  be  considered. 
Figure  6.2  shows  the  ray  paths;  Figure  6.2(a)  represents  paths  with  reflections 
and  Figure  6.2(b)  represents  equivalent  paths  by  the  method  of  images.  Four 
separate  paths  must  be  taken  into  account,  and  this  is.  sometimes  referred  to 
as  the  four-ray  problem.  A  program  has  been  written  to  solve  this  problem 
when  the  dimensions  are  specified  and  the  reflection  coefficients  given  for 
the  plane  on  each  side  of  the  mask.  A  complete  description  of  this  program, 
with  a  program  listing,  is  given  in  Appendix  A. 

We  have  calculated  the  pattern  propagation  factor  with  this  program  for 
a  representative  situation:  the  target  range  IS  km,  the  mask  10  km  from  the 
radar,  and  both  the  radar  and  the  edge  of  the  mask  are  100  m  above  the  re¬ 
flecting  plane.  For  each  radar  frequency,  four  cases  are  worked  out;  the 
one-way  propagation  loss  relative  to  free  space  (20  log  F)  is  calculated  for 
a  range  of  target  heights  0  to  250  m,  (1)  with  p  ■  0  on  both  sides  of  the 
mask,. (2)  fbr  reflection  on  the  left  side  only,  (3)  for  reflections  on  the 
right  side  only,  and  (4)  for  reflections  on  both  siies.  These  four  cases  are 
shown  for  an  L-band  radar  frequency  (1,3  GHz)  in  Figure  6.3.  When  p(left) 
and  p(right)  are  both  zero  we  get  the  familiar  Fresnel  diffraction  curve. 

When  reflection  occurs  on  the  left  only,  we  get  a  hybrid  situation  with  lobe 
structure  above  200  m  where  the  direct  and  reflected  waves  interfere,  and 
diffraction  of  both  waves  into  the  shadow  zone  with  little  change  in  their 
relative  phase.  When  reflection  occurs  on  the  right  only,  we  find  interfer¬ 
ence  lobing  between  the  direct  and  reflected  wave  with  increasingly  deep 
nulls  as  the  target  moves  into  the  shadow  zone.  We  can  think  of  the  knife- 
edge  as  the  source  of  the  diffracted  wave  in  the  shadow  zone,  and  the  lobe 
width  is  determined  by  the  height  of  the  knife-edge  above  the  reflecting 
plane.  When  reflection  occurs  on  both  sides  of  the  mask,  the  propagation 
combines  the  characteristics  of  the  two  cases  of  single  reflection. 

A  similar  computation  made  for  a  VHF  frequency  (170  MHz)  is  shown  in 
Figure  6.4;  the  geometrical  arrangement  is  identical  to  that  of  Figure  6.3. 
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Fig.  6.2  Diagram  of  the  four-ray  problem,  (a)  shows  the  direct  and  reflected  rays 
and  (b)  the  equivalent  construction  with  images. 


Fig.  6.3  Propagation  over  a  knife-edge  on  a  plane  at  a  frequency  of  1.3  GHz. 
The  propagation  geometry  is  shown  in  the  bottom  diagram;  the  four  graphs  show 
the  cases  of  perfect  reflection  or  no  reflection  on  each  side  of  the  knife- 


The  target-height  scale  for  this  radar  frequency  goes  only  up  to  250  m,  not 
high  enough  to  show  the  lobe  pattern  which  is  produced  by  reflections  in  the 
foreground  of  the  radar.  With  reflections  only  on  the  left  side  of  the  mask, 
the  signal  strength  is  increased  for  heights  from  150  to  250  m  and  decreased 
for  heights  below  100  m  as  compared  to  the  diffraction-only  case.  With 
reflections  on  the  target  side  of  the  mask,  one  finds  lobes  in  the  shadow 
zone  which  result  from  interference  between  waves  from  the  virtual  source  at 
the  knife-edge  as  in  the  L-band  case.  The  dotted  curves  in  Figure  6.4  repre¬ 
sent  the  diffraction-only  caje  for  comparison. 

6.4  Propagation  Over  a  Spherical  Earth 

The  problem  of  radio  propagation  over  a  smooth  spherical  earth  is  discussed 
in  detail  by  KERR  (1951),  by  Domb  and  Price  (1946),  and  by  many  others.  We 
describe  here  the  solution  of  this  classical  problem;  the  references  cited 
above  should  be  consulted  for  a  complete  discussion  of  the  problem.  The 
radius  of  the  sphere  is,  of  course,  taken  to  be  the  effective  earth  radius, 
i.e.,  radius  of  curvature  of  the  terrain  profile  corrected  for  atmospheric 
refraction.  For  a  radar  antenna  at  height  the  range  r^  to  the  horizon  is 
given  by  and  looking  from  the  target  toward  the  radar,  the  range 

r_2  to  the  horizon  is  ft"e,  where  hg  is  the  target  height.  Hence,  the 

range  r  at  which  the  line-of-sight  from  radar  to  target  is  just  tangent  to 
the  earth  (i.o.,  the  lowest  line-of-sight)  is  given  by: 

r.  +  r.  ■  /nrrrx  +  /2ettt 

i  i  i  e  2  e 

These  results  may  be  derived  from  simple  geometrical  considerations. 

The  solution  to  the  problem  will  be  considered  in  three  separate  regions. 
Figure  6.5  shows  how  these  regions  are  defined  in  terms  of  the  lowest  line- 
of-sight.  In  the  interference  region  the  propagation  is  determined  by  the 
interference  of  the  direct  and  reflected  waves.  The  path  difference  is 
computed  from  geometrical  considerations  and  the  amplitude  of  the  reflected 
wave  is  corrected  for  the  increased  divergence  of  the  rays  produced  by 
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INTERMEDIATE  REGION 


Pig.  6.S  Propagation  region*  over  a  smooth  spherical  earth. 


50 


reflection  from  a  convex  surface.  The  boundary  of  the  interference  region 
may  be  taken  as  the  line  along  which  the  path  difference  for  the  direct  and 
reflected  waves  equals  X/8.  The  diffraction  region  includes  points  that  lie 
sufficiently  far  below  the  tangent  ray  for  the  field  to  be  represented  by  a 
single  mode;  in  the  intermediate  region  all  modes  (the  complete  residual - 
series  solution)  are  required  to  describe  the  field.  Domb  and  Price  (1846) 
give  a  convenient  method  for  obtaining  the  pattern  propagation  factor  along 
the  tangent  ray  and  within  the  diffraction  region. 

We  have  programmed  the  interference-region  computations  of  the  pattern 
propagation  function  F;  a  description  of  this  program  is  contained  in  Appendix 
b.  As  an  example  of  this  calculation,  we  show  in  Figure  6.6  the  distributions 
of  values  of  at  various  ranges,  for  a  transmitting  antenna  at  a  height  of 
30  m  above  a  smooth  spherical  earth  (reflection  coefficient  ■  0.9  and  K  ■  4/3) 
calculated  for  wavelength  X  *  23  cm.  We  see  the  characteristic  lobe  structure 
represented  in  this  plot,  Note  that  the  lobes  and  peaks  line  up  along  a  straight 
line  in  Fig.  6.6.  The  lower  boundary  of  the  interference  region,  where  the 
path  difference  of  the  direct  and  reflected  wave  equals  A/8,  also  lies  along 
a  straight  line  in  Fig,  6.6,  The  figure  shows  values  of  F  which  vary  between 
1.9  (peaks)  and  0.1  (nulls).  The  pattern  of  the  transmitting  antenna  is  as¬ 
sumed  to  effectively  cut  off  the  signal  in  the  upper  left  corner. 

For  low-flying  aircraft  we  must  consider  propagation  loss  below  the 
interference  region.  Figure  6.7  shows  plots  of  F  vs.  altitude  for  a  target 
range  of  30  km  and  an  antenna  height  of  30  m.  The  lowest  line-of-sight  corre¬ 
sponds  to  a  target  height  of  3.25  m.  For  this  figure  we  assumed  p  ■  0.7  and 
K  ■  4/3;  the  wavelengths  represent  VHF  (X  ■  200  cm),  L-band  (X  ■  23  cm),  and 
X-band  (X  ■  3.3  cm).  The  boundary  of  the  interference  region  is  indicated  for 
each  wavelength.  The  boundary  is  at  the  following  target  heights:  170  tn  for 
X  ■  200  cm,  39  m  for  X  «  23  cm,  and  15  m  for  >3.3  cm.  For  the  range  of  al¬ 
titudes  in  Figure  6.7,  almost  the  entire  curve  for  VHF  wavelengths  lies  in  the 
intermediate  region,  while  at  X-band  most  of  the  curve  lies  in  the  interference 
region,  Note  that  the  depths  of  the  nulls  in  the  X-band  curve  decrease  as  the 
altitude  decreases.  This  is  a  consequence  of  the  divergence  factor  mentioned 
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Fig.  6.6  The  space  distribution  of  the  pattern  propagation  function  for  the 
interference  region  over  a  spherical  earth.  Parameter*  are  given  in  the  text 


Pig*  6.7  The  pattern  propagation  function  vs.  target  altitude  in  the  inter¬ 
ference  region  and  below.  Three  wavelengths  are  represented:  3.3,  23  and  200  cm; 
the  X/8  boundary  of  the  interference  region  is  indicated  for  each  wavelength  by 
the  change  from  solid  to  dotted  line.  The  lowest  lino-of-sight  corresponds  to 
a  target  height  of  3.2S  m. 
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earlier.  As  a  result  of  diffraction  by  a  smooth  spherical  earth,  we  expect 
that  the  VHF  wavelengths  will  suffer  a  much  greater  propagation  loss  than 
shorter  L-band  and  X-band  wavelengths  for  very  low-flying  targets.  These 
predictions,  however,  apply  to  propagation  over  a  highly  reflective  earth;  for 
example,  over  a  calm  sea.  We  know  of  no  complete  solution  for  the  interference 
region  that  specifically  includes  the  effect  of  surface  roughness,  but  in  the 
interference  region  one  can,  in  principle,  take  roughness  into  account  by 
assigning  a  reflection  coefficient.  In  the  spherical-earth  program  described 
in  Appendix  B,  the  reflection  coefficient  is  represented  by  a  parameter  that 
may  be  specified. 

7.  SUMMARY 

Low-altitude  propagation  at  frequencies  above  100  MHz  is  controlled  by  . 
atmospheric  refraction  and  by  diffraction  and  reflection  from  the  terrain  over 
which  the  waves  travel.  Where  these  phenomena  can  be  taken  into  account 
quantitatively,  we  can  predict  accurately  the  propagation  loss  as  a  function 
of  frequency  over  specific  terrain.  The  underlying  physical  principles  are 
understood,  but  a  number  of  difficulties  remain  before  the  low-altitude- 
propagation  problem  can  be  considered  solved, 

Refraction  plays  an  important  role  in  the  propagation  of  radar  waves  over 
large  bodies  of  water,  particularly  over  the  ocean  in  lower  altitudes,  Electro¬ 
magnetic  waves  may  be  trapped  in  surface  layers  of  air  containing  a  concentra¬ 
tion  of  water  vapor  and  an  inversion  of  the  temperature  gradient.  Greatly  ex¬ 
tended  radar  ranges  in  such  ducts  are  encountered  frequently  in  maritime  en¬ 
vironments.  On  the  other  hand,  over  continental  land  masses  ducting  is  in¬ 
frequent.  Nevertheless,  over  areas  of  desert  or  snow  cover,  nocturnal  tempera¬ 
ture  inversions  due  to  radiational  cooling  may  significantly  affect  propagation. 
But  usually  the  atmosphere  near  the  surface  over  land  has  a  refractivity  that 
decreases  nearly  linearly  with  increasing  altitude.  This  situation,  encountered 
roughtly  80  percent  of  the  time,  can  be  modeled  by  increasing  the  effective 
earth  radius  by  a  factor  which  varies  between  one  and  two,  a  typical  value 
being  4/3,  Taking  refraction  into  account  is  straightforward  in  this  case, 
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and  the  resulting  effect  on  radar  coverage  is  comparatively  small  and  in¬ 
dependent  of  radar  frequency. 

Diffraction  phenomena  on  the  other  hand  are  of  major  importance  in  low- 
altitude  propagation,  and  we  must  distinguish  between  two  cases,  In  the  first 
case  specific  irregularities  in  the  terrain  profile  produce  diffraction:  these 
terrain  features  can  be  represented  by  knife-edges,  or  in  some  cases  by  cylinders. 
In  the  second  case  diffraction  is  produced  by  the  curvature  of  the  smooth 
spherical  earth,  and  this  case  applies  to  propagation  over  oceans  or  very 
smooth  plains.  In  both  cases  the  propagation  is  strongly  dependent  on  radar 
frequency:  the  boundary  between  the  illuminated  region  and  the  shadow  becomes 
increasingly  sharp  as  the  frequency  increases.  At  lower  frequencies  diffraction 
by  hills  or  ridges  alters  the  field  strength  in  a  significant  region  above  and 
below  the  lowest  line  of  sight.  Diffraction  reduces  the  signal  above  the  mask 
and  at  the  same  time  directs  some  energy  into  the  shadow  region  below  the 
mask.  Consequently,  at  lower  frequencies  several  diffracting  features  in  the 
terrhin  profile  will  usually  contribute  to  the  propagation  loss,  This  presents 
a  problem  that  has  yet  to  be  completely  solved:  we  lack  an  adequate  model  for 
more  than  two  knife-edges  in  series.  Furthermore,  we  do  not  know  the  extent 
to  which  the  rigorous  solution  to  the  problem  of  propagation  over  a  smooth 
spherical  earth  is  applicable  when  the  terrain  is  not  perfectly  smooth.  These 
problems  require  thorough  investigation.  For  example,  we  should  be  able  to 
specify  the  conditions  under  which  the  multiple-diffraction  case  goes  over  to 
the  spherical  earth  case  as  the  terrain  roughness  decreases. 

The  effect  of  coherent  reflection  from  terrain  must  be  taken  into  account 
for  propagation  paths  over  water  and  smooth  plains,  particularly  at  lower  fre¬ 
quencies.  The  lobe  structure  produced  by  reflections  from  a  plane  gives  a  null 
at  the  horizon,  and  the  depth  of  this  null  depends  on  the  reflection  coefficient. 
So  coherent  reflection  from  terrain  can  alter  low-altitude  propagation  signifi¬ 
cantly  depending  on  the  value  of  the  reflection  coefficient.  If  the  terrain 
is  perfectly  smooth,  the  conductivity  and  dielectric  constant  of  the  surface 
determine  the  reflection  coefficient,  but  for  real  terrain  the  roughness  must 


b«  taken  into  account.  For  barren  terrain  the  Gaussian  terrain  model  may  be 
used  to  estimate  the  effect  of  terrain  irregularity,  but  this  model  has  not 
been  verified  for  very  low  grating  angles.  Furthermore,  the  reflection  prop¬ 
erties  of  terrain  with  vegetative  cover  are  almost  unknown  for  grazing  angles 
below  one  degree.  Thus,  extensive  measurements  over  barren  and  vegetated 
terrain  must  be  made  to  develop  a  valid  model  for  reflection  effects. 

The  standard  models  of  propagation  described  here  involve:  (1)  propaga¬ 
tion  over  a  plane  with  arbitrary  reflection  coefficient,  (2)  propagation  over 
a  sphere  in  a  region  above  the  horizon  where  the  direct  and  reflected  waves 
interfere,  (3)  propagation  over  a  knife-edge  on  a  reflecting  plane,  and  (4) 
propagation  over  successive  knife-edges  treated  approximately  by  the  method 
of  Deygout  (1966).  In  many  cases  these  models  are  adequate  for  predicting 
propagation  loss  provided  the  reflection  coefficients  are  known. 


APPENDIX  A 


FOUR-RAY  PROPAGATION  MODEL 

This  model  calculates  the  pattern  propagation  factor  F  for  radio  propa¬ 
gation  over  flat  terrain  on  which  a  knife-edge  obstruction  lies  perpendicular 
to  the  direction  of  propagation.  The  geometrical  arrangement  is  shown  in 
Fig.  A*  <  program  computes  F  for  a  sequence  of  values  of  Z 2,  the  height  of 
R  above 


Input 


Output 


Geometrical  distances  in  meters:  Z^,  dj, 
d2,  f,  and  the  initial  value,  final 
value,  and  step  for  Zg. 


Table  with  the  following  columns: 
Height  of  R;  Z2  (in  meters),  F 
(pattern  propagation  factor), 
and  20  log  F. 


Wavelength  A  in  meters 

Polarization:  Horizontal  (H)  or 
Vertical  (V) 

Ground  properties  on  left  and  right 

side  of  the  mask:  Relative  dielectric 
.constant  g ,  conductivity  o  (in  mhos/m) 


For  Rays  1  and  2 :  v ,  | E | ,  p 1 , 
♦j  and  v2,  |B2|,  p2,  <(>2 

For  Ray  4:  and  |fi^| 


(Note:  The  reflection  coefficient  on 
either  or  both  sides  may  be  taken 
as  zeio  if  e  -  a  -  888. ) 


In  case  the  surface  is  rough,  the  re¬ 
flection  coefficient  may  be  altered 
by  factors  REF  L  and  REF  R  on  the  left 
and  right  sides  of  the  mask. 


Reflections  on  the  left  and  right  sides  are  taken  into  account  by  including  the 
images  of  R  and  T  which  are  denoted  by  S  and  U,  respectively.  The  four  possible 
rays  to  be  considered  are  TR,  UR,  TS,  and  US,  which  are  designated  as  Rays  1,  2, 


3,  and  4,  respectively.  The  reflection  coefficients  are  designated  by  the 
subscripts  1  and  2  for  the  left  and  right  sides,  respectively.  Bach  reflection 
coefficient  T  is  given  by  a  complex  number  calculated  from  the  Fresnel  re¬ 
flection  formulas  in  a  subroutine,  FRBSNL,  and  Y  has  the  form 


where  the  magnitude  p  and  the  phase  lag  $  depend  on  the  grazing  angle,  the 
polarization,  and  the  electrical  constants,  e  and  o,  of  the  ground.  The  Fresnel 
integrals  C(v)  and  S(v),  are  evaluated  in  a  subroutine,  DCS,  where  v  ■  c 
with  c  equal  to  the  ray  clearance  over  the  knife-edge.  (When  the  ray  intersects 
the  mask,  c  is  taken  as  negative.)  The  pattern  propagation  factor  is  equal 
to  | F | ,  where 

4 

F  "  2  Bvexp  (MO 

k  ■  1  K 


and 


1  "  *1 


riA2 

r2A3 

rir2A4 


1 


2ir 


♦a  .  -  *2  +  T  CR1  ■  Ra> 


2ff 


rr 

tan 


*3  * 

*  H  ®4  •  Rl> 

l(C(v)  ♦  1/2) ^  +  (S(v)  +  1/2)2] 1/2 


it  +  tan 


1  /  S  +  1/2  N 

\ct ifr)  - 
1  /s  +  1/ 2  \ 

\rrm) 


tt/4  if  (C  +  1/2)  >  0 
-  tt/4  if  (C  +  1/2)  <  0 


This  computation  was  programmed  in  IBM  FORTRAN  IV  ANSI  by  Gerald  McCaffrey; 
a  copy  of  the  listings  for  this  program  as  well  as  a  sample  of  the  output 
follow  this  description.  The  sample  computation  appears  plotted  as  Fig.  6,4, 
and  all  the  input  parameters  are  listed  on  the  output. 
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non  cj  to  non  ^  non  t-1  non  non  r»oo 


PROGRAM  LISTINGS  >  FOUR -RAY  PROPAGATION  MODEL 
PROGRAMMED  BY  GERALD  MoCAFFREY 


VARIABLES  THAT  START  WITH  THH  LETTER  C.  ARB  COMPLEX 
COMPLEX  CT1 ,CT2 , CF,CE1 ,CE2 ,CB3 , CF4 , CTH , CTV 

REAL* 8  R1,R2,R3,R4,VC1,VC2,VC3,”C4,VS1,VS2,VS3,VS4,A,B,V1,V2,V3,V4 
REAL* 8  VTERM 
REAL* 4  LAMDA 
DIMENSION  POLAR (2) 

DATA  POLAR/ 'H  ','V  '/ 

PI-3.14159 
RAD-57.29578 

GET  SORT (2.0) 

SQRT2— SQRT (2.0) 

READ  GEOMETRIC  PARAMETERS 

WRITE (6,1) 

FORMAT ( '  INPUT  Z1,D1,D2,F,Z2S,Z2S,DZ2' ) 

READ ( 5  ,  * )  Zl,Dl,D2,F,Z2S,Z2B,DZ2 

READ  IN  DIALECTIC  CONSTANTS  AND  FRESNEL  FACTORS  (LEFT  AND  RIGHT) 

WRITE (6 ,4) 

FORMAT ('  INPUT  PSIL, SL, PSIR, SR,RBFL ,REFR' ) 

READ ( 5 , * )  EL , SL , 8R , SR , REFL , RBFR 

READ  IN  POLARIZATION 

WRITE (6 ,2) 

FORMAT ( '  INPUT  POLARIZATION..  H— HORIZONTAL, V— VERTICAL' ) 

READ (5,3)  POLR 
FORMAT (A4) 

IF ( POLR. EQ. POLAR (1) )  IPOL-1 
IF (POLR. EQ. POLAR (2))  IPOL-2 

INITIALIZE  Z2 

Z2-Z2S-DZ2 
C 

C  SET  WAVELENGTH 

C 
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on  o  n  o  o  o  o  o  o  h*  o  fi  n  on  o  i-*  m  <©  qoo 


LAMDA-.302 

WRITB  OUT  THE  INPUTS  AND  SET  UP  HEAD RE  FOR  THE  PRINTOUT 

WRITS (I , 10 5  Zl, Dl, D2 ,F,22Sf Z2B,DZ2 ,LAMDA,POLR,BL, SR, SL, SR 
WRITB (lr 9) 

FORMAT ( '  ' ) 

WRITB (6  »10)  Zl'DlrD2,F,Z2S,Z?.B,DZ2,LAMDA,POLR,HL,BR,Sti,SR 
l'RBFL'RSFR 

FORMAT ( 1  ANNT  HT 1 ,  F6 . 2 » '  D1',F6.0,’  D2',F6.0,'  F\F6.1,15X, 

1'  TARGET  HT  FROM' ,F5.0,'  TO',FE.O,'  BY* ,F5.0 WAVELENGTH' ,F6 . 3 , 
3'  POLARIZATION  \A4, '  BPSILON  L',F8.3,'  EPSILON  R\F8.3,'  SIG  L\ 
4F8.3,'  SIG  R' ,F8 . 3/ » '  REFL' , F1U . 2 , '  RBFR' rF10.2) 

WRITB (lr 11) 

FORMAT (2X, 'TARGET  HT' ,2X, '  F  ' ,1X' '20LOG(F) ' ,IX, '  V(l) »,1X, 

1'  MAG(El) '  ,2X, '  RHOd)',lX,'  PHKDMX,'  T(2)MX,'  MAG(B2)',3X, 
2'  V ( 3 ) ' » IX » 1  MAG (B3) ' , 1X» '  RH0(2)',1X,'  PHI(2)',3X,'  V(4)\1X, 

3'  MAG (E4) ' ) 

CALCULATE  HO 

HOl‘*LAMDA*Dl*D2 

H02-D1+D2 

HO-S0RT(HOl/HO2) 

LOOP  ON  Z2  FROM  Z2S  TO  Z2B  BY  DZ2 
Z2-Z2+DZ2 

IF (Z2 iGT, Z2E)  GO  TO  1000 

CALCULATE  Rl'R2,R3'AND  R4 

R1*DSQRT( (Z2*1.0D0-Zl)**2+ (D1+D2) **2) 

R2-DSQRT( (-22*1. ODO-Zl) **2+ (D1+D2) **2) 

R3-R2 

R4*R1 

CALCULAT  Vl,V2,V3,AND  V4 

VTKRM» ( (Z2*l .0D0-Z1) *Dl) / (D1+D2) 

V1-S0RT2/H0* (Zl+VTBRM-F) 

VTBRM** ( (22*1 , 0D0+Z1) *D1) / (D1+D2) 

V2-SQRT2/HO* (-Z1+VTHRM-F) 

VTRRM* ( (-21*1 . ODO-Z?) *D1) / (D1+D2) 

V3-8QRT2/HO*  (Zl*VTBRM-F) 

VTBRM- ( (Z1*1.0D0-Z2) *D1) t  (D1+D2) 

V4-SQRT2/HO* (-Z1+VTBRM-F) 

CALCULATE  THE  GRAZING  ANGLES 
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c 

C  LEFT  BRAZING  ANGLE 

C 

IF (V2 ,GE , 0 )  GO  TO  200 
C 

C  RAY  INTERSECTS  MASK 

C 

PSIL-DATAN<(Z1+F)/D1*1.0D0) 

GO  TO  201 
C 

C  RAY  CLEAR  OF  MASK 

C 

200  PSIL-DATAN((22+Z1)/(D1+D2*1.0D0)) 

C 

C  RIGHT  GRAZING  ANGLE 

C 

201  IFCV3.GB.0)  GO  TO  205 
C 

C  RAY  INTERSECTS  MASK 

C 

PSIR-DATAN( (Z2+F)/D2*1.0D0) 

GO  TO  206 
C 

C  RAY  CLEAR  OF  MASK 

C  ■ 

205  PSIR«DATAN( (Z2+Z1) / (Dl+D2*l. 0D0) ) 

C 

C  CALL  THE  COMPLEX  REFLECTION  COEFFICIENT  SUBROUTINE 

C 

C  FOR  THE  LEFT 

C 

206  CALL  FRBSNL ( EL , LAMDA f  SL , PS I L  r CTH , CTV ) 

C 

C  MULTIPLY  BY  FACTORS 

C 

CTH«CTH*RBFL 

CTV-CTV*RBFL 

C 

C  SET  CT1 

C 

IF ( IPOL . EQ . 1 )  CT1-CTH 
IF (IPOL.EQ . 2)  CTl-CTV 
C 

C  886  8H0WS  NO  REFLECTION  ON  THE  LEFT 

C 

IF (EL. BQ. 888. AND. SL. BO. 888)  CT1«CMPLX(0.0 f 0.0) 

C 

C  FOR  THE  RIGHT 

C 

CALL  FRBSNL ( HR , LAMDA, SR ,PSIR, CTH , CTV) 

C 

C  MULTIPLY  BY  FACTORS 


! 
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c 

CTH«CTH*REPR 

CTV-CTV*RBPR 

C 

C  SBT  CT2 

C 

■  IP (IPOL.BQ.l)  CT2-CTH 
IF(IPOL.BQ,2)  CT2-CTV 
C 

C  t»88  SHOWS  NO  RKFLBCTION  ON  THE  RIGHT 

C 

IP (BR. EQ. 986. AND. SR. BO. 688)  CT2«CMPLX(0.0,0.0) 

C 

C  OBT  THE  REAL.  AND  IMAGINARY  PARTS  POR  VBRT  AND  HORZ 

C 

P1R-RBAL (CT1 ) 

P1I“AIMAG(0T1) 

P2R-REAL  (CT2  ) 

P21-AIMAGCCT2) 

C 

C  GET  THB  PHASE  LAG  IN  DBGRBBS 

C 

THBTA1-0 

THSTA2-0 

IP (P1I .NB . 0 . AND. P1R.NB. 0)THSTA1-ATAN2 (P1I , P1R) * (-HAD) 

IF (P2I .N3 . 0 , AND. P2R.NB. 0) THBTA2-ATAN2 (P2I.P2R)* ( -RAD) 

C 

C  GBT  THB  REFLECTION  COBP 

C 

PR1-CA3S (CT1 ) 

PR2-CABS(CT2) 

C 

C  CALCULATE  COMPLEX  P  (CP) 

C  BY  PIRST  CALCULATING  THB  COMPLEX  PARTS  CE1 ,CB2 , CB3 , AND  CB4 

C 

C  CALCULATE  CB1  USING  VI 

C 

CALL  DCS (VC1 , VS1 , VI) 

A<"l/SORT2*DSORT  (  (VC1+0 . 5)  **2+  (VSl+0 . 5) **2) 

IP (VC1+0 . 5 .GB. 0)  B—DATAN ( (VSI+0 , 5) / (VC1+0 , 5) ) -PI/4 
IP (VC1+0 . 5 ,LT.  0)  B-PI+D'\TAN({VSl-t-0.5)/(VCl*0.5))-PI/4 
Ai-A*DSIN(B) 

Bl-A*DCOS(B) 

CBl-CMPLX(Bl.Al) 

C 

C  CALCULATE  CB2  USING  V2 

C 

CALL  DCS(VC2,VS2,V2) 

A~1/SQRT2*DSQRT ( (VC2+0 .5)**2*-(VS2  +  0,5)**2) 

IF (VC2+0 . 5 ,GE. 0)  B-DATAN ( (VS2+0, 5) / (VC2+0 . 5) ) -PI/ 4 
IF (VC2+0 . 5.LT. 0)  B-PI+DATAN ( (VS2+0 . 5) / (VC2+0 . 5) ) -PI/4 
B-B+ (2*PI) /LAMDA* (R2-R1) 


63 


A1«A*DSIN(B) 

B1-A*DC0S(B) 

CB2-CMPLX ( B1 , A1 ) *CT1 

C 

C  CALCULATE  CE3  USING  V3 

C 

CALL  DCS(VC3,VS3,V3) 

A-l/SQRT2*DSQRT((VC3+0,5)**2<(VS3+0.5)**2) 

IF(VC3+0.5.GB.0)  B-nATAN((VS3+0.5)/(VC3+0.5))-PI/4 
IF (VC3+0 .5 .LT. 0)  B»PI+DATAN( (VS3+0.5) / (VC3+0.5) )-PI/4 
S-B+ (2*PI) /LAMDA* (R3-R1) 

A1-A*DSIN(B) 

Bl-A*DCOS(B) 

CE3-CMPLX(B1,A1)*CT2 

C 

C  CALCULATE  CS4  USING  V4 

C 

CALL  DCS(VC4,VS4,V4) 

A-l  /  SQRT2*D8QRT  ( ( VC4 + 0 . 5 ) *  * 2+ ( VS4 + 0 . 5 ) *  *  2 ) 

IF(VC4+0 . 5 .OB. 0)  B“DATAN( (VS4+0.5) / (VC4+0.5) )-PI/4 
IF (VC4+0 . 5 .LT. 0)  B-PI+DATAN( (VS4+0.5) / (VC4+0.5) J-PI/4 
Al«A*D8IKt  (B) 

Bl-A*DCOS (B) 

CB4-CT1*CT2*CMPLX(B1,A1) 

C 

C  '  NOW  ADD  TO  OBT  CF 

C 

CF-CB1+CB2+CB3+CH4 

C 

C  WRITS  OUT  ANSWERS 

C 

FO-CABS (CF) 

FLOG«20*ALOG10(PO) 

KlMAO-CABS(CBl) 

B2MAG-CABS (CB2) 

B3MAG-CABF (CB3) 

B4MAO-CABB(CB4) 

WRITE (1,99)  Z  2 , FO , FLOG , VI , BIMAG , PR1 , THBTA1 , V2 , B2MA0 , V3 , B3MAG , 
1PR2 ,THBTA2 , V4 , B4MAG 
99  FORMAT (PS .1 , 4 (2F7. 2,F10 ,3) .2F8.2) 

GO  TO  100 
1000  RETURN 
BND 


non  non 


SUBROUTINE  FOR  COMPLEX  REFLECTION  COEFFICIENTS 


SUBROUT INB  FRBSNL (Bl.WAVB, CONDUC , ANO , CTH , CTV ) 
COMPLEX  CAK , CTV , CTV1 , CTV 2 , CTH , CTHl , CTH 2 
C 

C  B1 . THB  DIBLBCTIC  CONSTATNT  (PROM  0  TO  100) 

C  LAMDA. . . .THB  WAVBLBNOTH  1H  METERS 

C  CONDUC... THB  CONDUCTIVITY  IN  MHOS/MBTBR 

C  ANO . THB  ANOLB  IN  RADIANS 

C 

C 

C  CALCULATE  THB  COMPLEX  CONSTANT 

C 

AKI— 60*WAVS*CONDUC 
CAK-CMPLX(Bl.AKI) 

CALCULATE  THB  VERTICAL  POLARIZATION 


CTV1-CAK*  SIN (ANO ) -CSQRT (CAK-COS ( ANO ) *  *  2 ) 
CTV2-CAK* SlN (ANO ) + CSQRT (CAK-COS ( ANO ) *  *  2) 
CTV-CTV1/CTV2 


CALCULATE  THB  HORIZONTAL 

CTHl-SIN(ANO) -CSQRT (CAK-COS (ANO) *»2) 
CTH2-S IN (ANO ) +CSQRT ( CAK-COS (ANO ) * *  2 ) 
CTH-CTH1/CTH2 
1000  RETURN 
BND 


SUBROUTINE  TO  EVALUATE  THE  FRESNEL  INTEGRALS 

SUBROUTINB  DCS  (C,S,X) 

IMPLICIT  REAL *6  (A-H,0-Z) 

DIMENSION  CC (12) ,DD(12) ,AA(12) ,BB(12) 

U-X 

PIH2-1.5707963268DO 

X-PIB2*X*X 

Z-DABS(X) 

IP(Z.NB.  0.D0)  OO  TO  1000 

C-O.DO 

S-O.DO 

x-u 

RETURN 

1000  CONTINUB 

IF(Z-4.0D0)3r3,4 
3  C-DCOS(Z) 

S-DSIN(Z) 

Z-Z/4.0D0 

DZ-DSQRT(Z) 
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AA(1)-0.159576914D+01 
AA(2)—0. 17020-05 
AA(3)— 0.68085688540+01 
AA(4)— 0.576361D-03 

AA(5)"0. 69206919020+01 

AA(6)  —  0.16898657D-01 

AA(7)"-0.305048566D+01 
AA(8)— 0.75752419D-01 
AA (9 5-0.85066378100 
AA(10)«-0 • 256390410-01 
AA(ll)— 0.15023096D0 
AA(12)"0 . 34404779D-01 
BB(1)"-0.33D-07 
BB ( 2) “0 .42553875240+01 
BB(3) n-0 . 9281D-04 
BB(4)--0. 778002040+01 
BB(5)— 0.9520895D-02 
BB( 6 )“0 . 5075161298D+01 
BB(7)— 0.138341947D0 
BB(8)"-0 • 13C3729124D+01 
BB(9)— 0.40334927600 

BB(10)“0 .7O2222016D0 

BB(ll)— 0.21619592900 

BB(12)"O.19547O31D-01 

ASUM-AA(l) 

BSUM-BB(l) 

DO  40  J-2,12 
ABUM-ASUM+AA(J)  +z**  (J-i) 
BSUM-BSUM+BB(J)*Z** (j-l) 
40  CONTINUE 

PC-DZ* (S*BBUM+C+A8UM) 
PS-DZ* (-C+BBUM+8+A8UM) 
C-PC 
S-PS 
00  TO  5 
4  D-DCOB(Z) 

S-DSXN(Z) 

Z"4 . ODO/Z 
CC(1)"0.0D0 
CC(2)— 0.249339570-01 
CC (3)"0 ,39360-05 
CC (4 )"0 . 5770956D-02 
CC(5)-0. 689B92O-03 
CC (6)“-0 . 9497136D-02 
CC (7)-0 . 11948809D-01 
CC (8 )»-0. 67488730-02 
CC (9)"0 . 24642D-03 
CC  (10)“0 . 2102967D-02 
CC(ll)— 0.121793D-02 
CC(12)-0. 2339390-03 
DD(1)-0.19947114D0 
DD (2) -0.230-07 
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DD<3)  — 0.9351341D-02 
DO (4) *0.230060-04 
DD(5)-0. 485146D-02 
DD(6)-0. 19032180-02 
DD<7)  — 0.17122914D-01 
DD(8)»0. 290640670-01 
DD(9J— 0.2792895BD-01 
DD(10)-0.16497308D-01 
00(11)— 0.6598515D-02 
DD(12)-0.838386D-03 
DSUM-DD(l) 

CSUM*CC  ( 1 ) 

DO  30  J-2,12 
CSUM-CSUM+CC ( J ) *  Z*  * ( J-l ) 
DSUM-DSUM+DD(J) *Z** (J-l) 
Z-DSQRT(Z) 

C-0 , 5D0+Z* (D*CSUM+S*DSUM) 
8-0. 5D0+Z* (S*C8UM-D*DSUM) 

x-u 

IF(U.OT.O.DO)  00  TO  6 
C— C 
8—8 
RBTURN 
END 
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APPENDIX  B 


THE  SMOOTH  SPHERICAL-EARTH  MODEL  FOR  THE  INTERFERENCE  REGION 


This  model  implements  in  IBM  FORTRAN  IV  ANSI  the  computation  of  propa¬ 
gation  in  the  interference  region  over  a  smooth  earth  as  described  in  KERR 
(1951) ,  pp.  112-122. 

For  a  specified  radar  wavelength,  transmitting  antenna  height,  and  earth 
reflection  coefficient,  the  program  described  here  computes  the  pattern  propa¬ 
gation  factor  for  a  sequence  of  target  heights  at  a  specified  target  range. 
Refraction  is  taken  into  account  by  specifying  the  earth-radius  factor  K;  the 
true  radius  is  taken  as  6373  km.  All  parameters  and  variables  are  desig¬ 
nated  by  the  symbols  used  in  that  reference  unless  specified  otherwise;  the 
input  and  output  quantities  are  as  follows: 

Output 

Table  with  the  following  columns: 

Target  Height 

Pattern  Propagation  Factor  F 

F 

20  log  F 
40  log  F 
Grazing  angle 
Distances 

R^  transmitter  to  specular  point 
R2  target  to  specular  point 

Divergence  factor 

AR  path  difference  between 
direct  and  reflected  waves 

First-Fresnel-Zone  size 

S-MAJ  semimajor  axis 

S-MIN  semiminor  axis  of  the 
first  Fresnel  zone 


Input 

Target  range  (km)  RNG 
Transmitter  antenna  height  (m)  Z1 
Earth  radius  factor  AK 
Wavelength  (cm)  WAVE 
Reflection  coefficient 
Magnitude  RHO 
Phase  PHAZ 
Target  height  (m) 

Starting  height  Z2S 
Ending  height  Z2E 
Incremental  step  DZ2 
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Note  that  the  output  contains  a  number  of  geometrical  quantities  that  axe 
useful  in  evaluating  the  model  predictions.  For  example,  knowing  the  path 
difference  AR,  one  can  locate  the  boundary  of  the  interference  region, 

AR  *  X/8.  The  dimensions  of  the  first  Fresnel  zone  are  calculated  approxi¬ 
mately  by  constructing  a  plane  tangent  to  the  earth  at  the  specular  point  and 
calculating  the  dimensions  of  the  semimajor  and  semiminor  axes  of  the  Fresnel 
ellipse  on  this  tangent  plane  (see  KERR,  1951,  p.  415,  Eqs.  (31)-(33)). 

This  computation  was  programmed  by  Gerald  McCaffrey;  a  copy  of  the  pro¬ 
gram  listings  and  a  sample  output,  which  includes  a  list  of  input  parameters, 
follow  this  description.  The  sample  output  is  plotted  for  X  ■  3.3  cm  in 
Fig.  6.7. 


PROGRAM  LISTINGS:  SPHERICAL -EARTH  MODEL, 
INTERFERENCE  REGION 


PROORAMMED  BY  OERALD  MoCAFFREY 


DATA  PI/3. 14159/ 

DATA  RAD/57.2957795/ 

DATA  RADIUS/6373. 0/ 

C 

C  RBAD  INPUTS ,  PROM  INPUT  PORM 
0 

WRITE (6 ,10) 

10  FORMAT ('  INPUT  TRAN  HT, WAVELENGTH, BARTH  RAD  FACTOR, RHO  AND  PHASE1) 
RBAD ( 5 , * )  Z l , WAVE , AK , RHO , PHAZ 
PZ1-Z1 

C 

C  CONVERT  PHAZ  TO  RADIANS 

C 

RPHZ-PHAZ/RAD 

C 

C  '  CALCULATE  AB  -  THE  BFPBCTIVB  BARTH  RADIUS 
C 

AB-RADIUS*AK 

C 

C  GET  RANGE  AND  TARGET  LOOPS  SETUP 

C 

WRITE (6, 204) 

204  PORMAT ( *  INPUT  RANGE  LOOP  I  RNGS , RNGE , DRNG 1 ) 

RBAD (5 , * )  RNGS, RNGE, DRNO 
WRITB(6 , 203) 

203  PORMATC  INPUT  TARGET  HEIGHT  LOOP  »Z2ST,Z2B,DZ2' ) 

RBAD (5,*)  Z2ST , Z2B,DZ2 
Z2S-Z2ST-DZ2 
C 

C  LOOP  ON  RANGE 

C 

RNG-ANGS-DRNG 
200  RNGmRNG + DRNG 

IP(RNG,OT.RNGB)  GO  TO  1000 
IC-0 
IC1-0 
C 

C  LOOP  ON  TARGET  HEIGHT 

C 

Z2S-Z28T-DZ2 
250  228-Z2S+DZ2 

IP(Z2S,OT,Z2B)  GO  TO  200 

Z2-Z28 

PZ2-Z2 

71 


QOO  to  o  noo  o  o 


CHECK  THAT  Zl.OB.22 
Z1-PZ1 

IF(Z1.GE.Z2)  00  TO  20 

IF  HOT  SWITCH  21  AND  22 

SAVE-21 
21*22 
Z2-SAVB 

IC-IC+1 
IC1-IC1+1 

CALCULATE  R2  -  DISTANCE  TARGET  TO  REFLECTION  (KM) 

A-AB* ( (Zl+Z2)/1000,0) 

B-(RNO/2.0)**2 
P-2. 0 /SORT (3  .0)  *SQRT  (A+B) 

ANG-ARC0S ( (Z2-Z1) *2*AH*RNO/ (1000.0*P**3) ) 

R2-RNO/ 2 . 0 ♦ P*COS ( ( ANO+PI ) / 3 . 0 ) 

R1-RNO-R2 
PR1-R1 
PR2-R2 

IF (PZ1 .0B. PZ2)  00  TO  25 
SAVE-PR1 
PR1-PR2 
PR2-SAVB 
25  Sl-Rl/SQRT(2*AB»21/1000,0) 

S2-R2/ SORT (2*AB* 22/1000 . 0 ) 

S-RN3/ (SORT(2*AB*Zl/1000.0)+BQRT(2*AB*Z2/1000.0) ) 

IP(Sl.LT.l.AND.B2.LT.l.AND.S,LT.l)  00  TO  30 

F-0 

F20-0 

F40-0 

X-0 

Rl-0 

R2-0 

00  TO  39 

30  T-SQRT (22/21) 

D1-4*S2*T*S1**2 
D2-S* (1-S1**2) * (1+T) 

D3-1+D1/D2 
D-l/SQRT (D3) 

Xl“  ( (1-S1**2)  +T*'*2*  (1-82** 2) )  /  (1*T**2) 

X2- (Zl+Z2)/(1000 , 0*RNO) 

X3«X1*X2 
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X=ATAM(X3) 

XDEG-X*RAD 
DR1-2*Z1*Z2/ (10*RNG) 

DR2-(1-Sl**2)*(l-S2**2) 

DR-DR1*DR2 

Fl«l+RHO**2*D**2 

F2-2 *RHO*D*COS  ( (2*PI*DR/WAVB)  4-RPHZ) 

F-SQRT(F1+F2) 

F20-20*ALOG10(F) 

F40»40*ALOG10(F) 

Z1P-Z1-R1**2/(2*AB) 

Z2P-Z2-R2**2/(2*AE) 

CG-1+2*DR/WAVE 

CH-1  + (Z1P+Z2P) **2/ (10*WAVB*RNG) 

CI-SORT(10*WAVB*RNG) 

CJ-1+2*Z1P* (Z1P+Z2P) / (10*WAVB*RNG) 

X0-RNG/2*CJ/CH 
IF ( PZ2 . GT . PZl )  XO-RNG-XO 
XC-RNG / 2  *  SORT ( CG ) /CH 
YC-C I / 2  * SORT ( CG /CH ) 

39  IF(IC.OT.l)  GO  TO  40 
WRITB(8,0999) 

8999  FORMAT(50X,'PROPl (12/11/78  VERSION)') 

WRITS (8 , 9000)  RNG , WAVE 

9000  FORMAT  Cl* ,5X, 'RANGE-  ',F4.1,‘  KM' ,70X, 'WAVELENGTH-  ' , F5 . 1 , '  CM') 
WRITB (6 , 9001)  PZl 

9001  FORMAT (5X, 'TRANSMITTER  HEIGHT-  ' ,F5.1,'  M' , 57X, 'RBPLBCTION  COEFFIC 
1IENT' ) 

WRITB (8 i 9002)  AK,RHO 

9002  FORMAT (5X, 'EARTH  RADIUS  FACTOR"  ' , F5. 2,58X, 'MAGNITUDE-  'rF5.3) 
WRITB (8 , 9003)  PHAZ 

9003  FORMAT ( 9  9X , ' PHASE-  ' ,F6.1,/) 

WRITE(8 , 9004) 

9004  FORMAT (2X, 'TARGET  HEIGHT ', 4X ,* PATTERN  PROPAGATION  FACTOR ', 3X ,' GRAZ 
1 ING ' , 0X , ' Rl ' , 7X , ' R2 ' , 7X , ' DI V ' , 5X , ' DBLTA-R ' , 6X , ' FIRST  FRESNEL  ZONE 
2') 

VJRITB(8  ,  9005) 

9005  FORMAT ( 7X , ' (M) ' ,10X, 'F' ,  5X,  ' 20 LOG (F) ' ,  2X, ' 40LOG (F) ' ,2X* ANGLE  (DBG) 
1' ,5X, ' (KM) ' ,5X,' (KM) ', 5X, ' FACTOR' , 5X , ' (CM) ' ,2X, ' CENTER (KM) ' , IX, 

2 ' S-MAJ (KM) ' ,  2X , ' S -MINOR (M) ' ) 

40  WRITE (6 , 9006)  PZ2 , F,F20 ,F40 ,XDBG , PR1 , PR2 ,D, DR,XO,XC ,YC 

9006  FORMAT (F10 . 2 ,F13 . 2 ,F10 . 2 ,F10 , 2 , F10 . 2 , F13 .2,F9.2,F11.3,F9,2,F9.3,F9 
1.3, FI 3. 2) 

WRITB (1 ,45)  RNG,PZ2,F20,F40 
45  FORMAT (4F12. 4) 

IF(IC1. BQ. 4)  WRITB (8 , 9007) 

9007  FORMAT ('  ') 

IF (ICl , BQ. 4)  IC1-0 
GO  TO  250 
1000  RETURN 
END 
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SAMPLE  OUTPUT:  SPHERICAL  -  EARTH  MODEL,  INTERFERENCE  REGION 
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APPENDIX  C 


BIBLIOGRAPHY:  RADAR  PROPAGATION  AT  LOW  ALTITUDES 

''This  bibliography  covers  the  subject  of  radio  and  microwave  propagation, 
VHP  through  X-band,  near  the  surface  of  the  earth.  References  are  indexed 
under  the  headings  of  Books,  Journal  Articles,  and  Technical  Reports  and  cross 
referenced  in  a  Sub j  ect  Index .  Although  many  of  the  references  cited  here 
deal  with  one-way  propagation,  the  results  are  directly  applicable  to  the 
radar  (two-way)  case  as  explained  in  Sec.  2.  We  assembled  this  bibliography 
as  we  surveyed  the  literature  on  the  effects  of  propagation  on  ground-radar 
performance  against  low-flying  aircraft.  All  listed  references  were  examined 
and  found  relevant  to  this  general  subject;  the  references  judged  most  useful 
are  marked  with  an  asterisk.  In  striving  for  completeness  we  have  included 
some  references  of  marginal  value,  but  all  are  worthy  of  examination  if  their 
title  or  their  listing  in  the  Subject  Index  suggests  that  they  may  be  of  inter 
est.  References  in  the  text  of  this  report  and  in  the  Subject  Index  are 
labeled  by  the  last  names  of  authors  and  the  year  of  publication.  The  Author 
Index  contains  the  complete  reference  information.  Major  headings  in  the 
Subject  Index  are  listed  alphabetically  below. 
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2 .  Frequency  Bands 
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4.  Microwave  Links 

5.  Mobile  Communication 
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7.  Ocean,  Propagation  Over 

8.  Reflection  Effects 
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11.  Spherical  Earth,  Propagation  Over 
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